Preemptive nematic order, pseudogap, and orbital order in the iron pnictides 
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Starting from a microscopic itinerant model, we derive and analyze the effective low-energy model 
for collective magnetic excitations in the iron pnictides. We show that the stripe magnetic order 
is generally preempted by an Ising-nematic order which breaks Ca lattice symmetry but preserves 
0(3) spin-rotational symmetry. This leads to a rich phase diagram as function of doping, pressure, 
and elastic moduli, displaying split magnetic and nematic tri-critical points. The nematic transition 
may instantly bring the system to the verge of a magnetic transition, or it may occur first, being 
followed by a magnetic transition at a lower temperature. In the latter case, the preemptive nematic 
transition is accompanied by either a jump or a rapid increase of the magnetic correlation length, 
triggering a pseudogap behavior associated with magnetic precursors. Furthermore, due to the 
distinct orbital character of each Fermi pocket, the nematic transition also induces orbital order. 
We compare our results to various experiments, showing that they correctly address the changes in 
the character of the magneto-structural transition across the phase diagrams of different compounds, 
as well as the relationship between the orthorhombic and magnetic order parameters. 



I. INTRODUCTION 

The current interest in iron-based superconductors 
(FeSCs) lies not only in the superconducting pairing 
mechanism, but also in the peculiar normal-state prop- 
erties of these materials (for reviews, see [1]). Of par- 
ticular interest is the fact that in weakly doped iron 
pnictides the stripe spin-density wave order at T < T/v, 
with ordering vectors (0, it) or (tt, 0) in the 1-Fe Bril- 
louin zone [2, 3], is often preceded by a "nematic" phase 
with broken C4 tetragonal symmetry but unbroken 0(3) 
rotational symmetry. The emergence of such a phase 
is not only manifested by a tetragonal to orthorhombic 
transition at T s > TV, but also by the onset of signifi- 
cant anisotropics in several quantities [4], such as dc re- 
sistivity [5, 6], optical conductivity [7, 8], local density 
of states [9], orbital occupancy [10], uniform suscepti- 
bility [11], and the vortex core in the mixed supercon- 
ducting state [12] . The fact that the spin-density-wave 
and structural transition lines follow each other across 
all the phase diagrams of 1111 materials (JTFeAsO, with 
rare-earth R) and 122 materials (AeFe2As2,with alkaline- 
earth Ae), even inside the superconducting dome [13, 14], 
prompted researchers to propose that these phases are in- 
timately connected. The interplay between magnetic and 
structural transitions in FeSCs is also quite rich: while 
in 1111 materials the two transitions are second-order 
and split (T s > T/v), in most of the 122 materials they 
seem to occur simultaneously or near-simultaneously at 
small dopings, but clearly split above some critical dop- 
ing - x ~ 0.022 in Ba (Fei_ x Co2;) 2 AS2, see [15, 16], and 
x sa 0.039 in Ca (Yei^ x Go x ) 2 AS2, see [17]. 

Early theoretical proposals explored two alternative 



scenarios for the nematic order. One scenario, which bor- 
rows concepts from the manganites, is that the driving 
force is orbital ordering - it induces the structural transi- 
tion and triggers the magnetic transition at a lower tem- 
perature by renormalizing the exchange constants [18- 
27]. An alternative, magnetic scenario, which borrows 
concepts from the studies of the J\ — J2 model of lo- 
calized spins [28], is build upon the fact that the or- 
der parameter manifold for the stripe magnetic order is 
0(3) x Z2, with the Z2 Ising symmetry differentiating be- 
tween the two possible ordering vectors (0, it) and (it, 0). 
This scenario explores the possibility that the Z2 Ising 
variable orders before the 0(3) rotational symmetry is 
broken, leading to an intermediate phase with a broken 
tetragonal symmetry but no long-range magnetic order. 
The Ising order makes the Fermi pockets at (0, tt) and 
(•7T, 0) non-equivalent and induces orbital order, since the 
two electron pockets are constituted of different orbitals 
[29, 30]. Furthermore, the breaking of the Z2 symme- 
try also makes the x and y directions inside the unit 
cell inequivalent, inducing a structural instability via the 
bi-linear coupling between the Z2 order parameter and 
the orthorhombic distortion (see below). For this rea- 
son, the state with broken Ising symmetry is often called 
an Ising-nematic phase. The magnetic scenario was first 
applied to localized or nearly-localized spin models for 
Fe-pnictides [31-33], but was later extended to itinerant 
systems in a phenomenological way [34, 35]. 

Since orbital order and Ising-nematic order break the 
same symmetry, one order generates the other, making 
the experimental distinction between the two scenarios 
rather subtle. The same is true if one wants to distin- 
guish between itinerant and localized magnetic models, 
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Figure 1: (Color online) Two characteristic phase diagrams 
of the system behavior as function of temperature and ne- 
matic coupling. Panel (a) is for strongly anisotropic, quasi- 
two dimensional systems, while panel (b) is for moderately 
anisotropic, less quasi-two dimensional systems. Red (light 
gray) and blue (dark gray) curves represent magnetic and 
Ising-nematic transitions, respectively. The arrows show how 
the inverse nematic coupling changes as function of doping, 
pressure, and lattice stiffness. Here and in the other figures, a 
solid (dashed) line denotes a second-order (first-order) transi- 
tion, and a double-dashed line indicates a simultaneous first- 
order transition. The two solid points mark the positions of 
the nematic and magnetic tri-critical points. The difference 
between the two phase diagrams is in the system behavior at 
intermediate nematic coupling. In (a), the system displays a 
first-order nematic transition followed by a second-order mag- 
netic transition at lower temperatures. In (b), the system first 
undergoes a second-order nematic transition, and then, at 
smaller temperatures, a meta-nematic transition which trig- 
gers a first-order magnetic transition. At the meta-nematic 
transition, the nematic order parameter jumps between two 
finite values. 



which we will address as well. In both situations, the only 
real way to distinguish between different approaches is by 
explicit calculations followed by comparison to the avail- 
able experimental results. In this regard, the richness 
of the phase diagrams of the iron pnictides is an impor- 
tant cornerstone, because one has to explain not only the 
very existence of the nematic phase, but also how this 



transition changes as a function of doping, pressure, and 
material. In addition, one has to address the feedback 
effect from the nematic order on fermions. In particu- 
lar, there is growing volume of experimental evidence of 
pseudogap-like electronic behavior in the iron pnictides, 
which in some cases seems to emerge at the same tem- 
perature at which nematic order sets in [36-39]. 

In this paper, we argue that the changes in the charac- 
ter of the nematic transition with doping, pressure, and 
alkaline-earth substitutions, as well as the development of 
the pseudogap and other experimentally detectable fea- 
tures, can be understood within a magnetic scenario for 
an itinerant fermionic model. We depart from a multi- 
band model of interacting fermions and derive (instead 
of assuming) the Ginzburg-Landau (GL) effective action 
for the two low-energy collective 0(3) magnetic degrees 
of freedom Ax and Ay associated with the ordering vec- 
tors (tt, 0) and (0, 7r), respectively. This action contains, 
besides other terms, the term j(A| — A y ) 2 which de- 
scribes fluctuations in the Ising-nematic channel. We ex- 
tend the original 0(3) model to N field components and 
study the limit N = oo within mean-field theory, and ar- 
bitrary N by the renormalization group (RG) technique. 
We find different types of system behavior depending on 
the strength of the dimensionless parameter a a j _1 
and on the degree of magnetic anisotropy. In all cases, 
however, we find two distinct multi-critical points in the 
phase diagram, namely a magnetic tri-critical point and 
a nematic tri-critical point. 

In the case of strongly anisotropic, quasi-2D systems, 
the nematic tri-critical point precedes the magnetic one, 
resulting in three types of behavior - see Fig. 1(a). (i) For 
small a (largest nematic coupling g) we find a strong first- 
order Ising-nematic transition at T s , accompanied by a 
discontinuous jump of the magnetic correlation length £ 
and a jump of the magnetic order parameter to a finite 
value. In this case, the stripe magnetic order emerges 
simultaneously with the Ising-nematic order, via a first- 
order magnetic transition (TV = T s ). but because £ 
jumps, there are no critical magnetic fluctuations above 
the transition, (ii) For intermediate a we find a first- 
order preemptive Ising-nematic transition accompanied 
by a discontinuous increase of the magnetic correlation 
length to a larger but still finite value, such that the stripe 
magnetic order does not develop at T s and emerges only 
at a smaller T. In this case, the magnetic transition is 
second-order, (iii) For large a (small nematic coupling 
g) we find a second-order preemptive Ising-nematic tran- 
sition followed by a second-order magnetic stripe transi- 
tion at a smaller T. Near the nematic tri-critical point 
we find that the Ising-nematic order emerges very rapidly 
such that, in practice, the second-order Ising transition 
is hard to distinguish from a first-order transition. The 
behavior of the nematic and magnetic order parameters 
for the phase diagram in Fig. 1(a) is shown in Fig. 9(d), 
(c), and (b) for regimes (i), (ii), and (iii), respectively. 

In the case of moderately anisotropic, less quasi-2D 
systems, the positions of the nematic and magnetic tri- 



critical points can be reversed, resulting in a new be- 
havior in the regime of intermediate a, see Fig. 1(b). 
The regimes (i) and (iii) are still present, when both 
the nematic and magnetic transitions are first-order and 
second-order, respectively. However, in regime (ii) the 
upper Ising-nematic transition is second-order and at a 
lower T there is a meta-nematic transition (i.e. the ne- 
matic order parameter undergoes a finite jump from one 
non-zero value to the other), which induces a first-order 
magnetic transition. The behavior of the nematic and 
magnetic order parameters for the phase diagram in Fig. 
1(b) is shown in Fig. 15(c), (b), and (a) for regimes (i), 
(ii), and (iii), respectively. 

As shown in Figs. 1, we find from our microscopic 
model that a decreases with pressure but increases with 
increasing lattice stiffness and electron doping, i.e. larger 
electron doping results in a larger value of a. Our results 
then predict that the Ising-nematic and magnetic transi- 
tions split upon electron doping, but tend to remain si- 
multaneous and first-order upon pressure, in agreement 
with the experimental data for FeSCs. We also show 
that the nematic order parameter couples linearly to the 
orbital polarization and lattice distortion, and hence ne- 
matic order generates orbital and structural order. We 
argue, however, that the sign of the orbital order may 
differ between hole and electron-doped materials. 

We also argue that, for the phase diagram in Fig. la, 
the discontinuous (or nearly-discontinuous) increase of 
the magnetic correlation length at T s in the regimes (ii) 
and (iii) greatly increases the strength of thermal mag- 
netic fluctuations. Once enhanced, these fluctuations ac- 
count for spectral weight redistribution (magnetic pre- 
cursors), what gives rise to pseudogap behavior in the 
fermionic spectral function and other observables. We 
present more detailed comparison with the data later in 
the paper and also compare our results with earlier stud- 
ies of Ising-nematic order. 

In the itinerant picture which we adopt here, the na- 
ture of the Ising-nematic phase has a clear interpretation 
in terms of magnetic fluctuations. The system has two 
degenerate stripe magnetic ground states with ordering 
vectors (n,0) and (0,7r), described by the two order pa- 
rameters Ax and Ay. At high temperatures (Ax) = 
(Ay) = 0, and the fluctuations of each order parame- 
ter have equal strength, i.e., (A x ) = (Ay). The Ising- 
nematic phase emerges when fluctuations associated with 
one of the ordering vectors become stronger than the 
other, (A 2 X ) ^ (Ay), while still (A x ) = (Ay) = 
(see Fig. 2). Because there are two possible choices, 
\A X ) > (Ay) or (A x ) < (Ay), the symmetry which 
breaks down at the Ising-nematic phase transition is a Zi 
symmetry. Once magnetic fluctuations around (jr, 0) be- 
come stronger (or weaker) than the fluctuations around 
(0, 7r), the equivalence between x and y directions inside 
the unit cell breaks down. The Ising-nematic transition 
then triggers the transition from a tetragonal to an or- 
thorhombic phase, and also imposes orbital order because 
the Fermi pockets centered at (tt, 0) and (0, tt) have differ- 
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Figure 2: (Color online) (upper panel) The band-structure 
with a circular hole pocket at V and two electron pockets at 
X and Y . The Brillouin zone contains one Fe atom, (lower 
panels) Static magnetic susceptibility Xq across the Brillouin 
zone for different temperatures. At high temperatures, fluc- 
tuations near the two stripe magnetic ordering vectors are 
equally strong, (A x ) = (Ay). Above the magnetic ordering 
temperature Tjv but below the Ising-nematic ordering tem- 
perature T s , fluctuations associated with one of the stripe 
states become stronger (in the figure, (A x ) > (Ay)) and the 
tetragonal symmetry is broken inside the unit cell. Stronger 
fluctuations around one ordering vector yield stronger inten- 
sity and narrower peaks. 



ent orbital character. Furthermore, this anisotropy in the 
spectrum of magnetic fluctuations leads to anisotropic 
scattering of electrons, resulting in different in-plane re- 
sistivities along x and y directions [40]. 

The structure of the paper is the following. In Section 
II we derive the effective action for the nematic order pa- 
rameter using a two-stage Hubbard-Stratonovich trans- 
formation. In Sec. Ill we extend the effective action to 
arbitrary N, where N is the number of components of the 
magnetic order parameter (N = 3 in the original 0(3)- 
isotropic model). We solve for the nematic order in the 
mean-field approximation, justified in the N — > oo limit, 
obtain different regimes of system behavior and discuss 
each regime in some detail. In Sec. IV we study the effec- 
tive action at a finite N using the renormalization group 
(RG) technique. We find that the system behavior for 
N = 3 is the same as for large N, but there is a change 
in the RG flow structure at N < 2. In Sec. V we discuss 
the interplay between Ising-nematic, orbital, and struc- 
tural order, and show how the nematic order gives rise 
to pseudogap behavior despite the fact that the Z 2 or- 
der parameter has zero momentum and does not by itself 
reconstructs low-energy fermionic states. In Sec. VI we 
compare our theory with the experiments and with other 
theoretical works on the nematic order. We present our 
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conclusions in Sec. VII. 



II. MICROSCOPIC MODEL AND THE 
EFFECTIVE ACTION 

We start from the minimal four-band model with two 
hole pockets T at the center of the Brillouin zone and two 
electron pockets X and Y at Qi = (n, 0) and Q2 = (0, n), 
respectively (see Fig. 2). We follow Ref. [41] and con- 
sider that only one of the two hole pockets participate in 
the SDW state. The extension to the case when both 
T pockets are involved in the SDW reconstruction is 
straightforward and does not change the analysis below 
except for the renormalization of the couplings. 

For simplicity, we consider parabolic dispersions with 

k 2 k 2 k 2 

£r,k = £0 - 2^ - £x,k+ Ql = -eo + 5^ + 2^ - /*> 

k 2 k 2 

and £y,k+Q 2 = so + ^ + 2^ ~ Mi where m, denotes 
the band masses, £0 is the offset energy, and n is the 
chemical potential. Near the Fermi energy and for small 
ellipticity, the dispersions can be approximated by £r,k = 
£x,k+Q! = £ - S + 5 2 cos29, £y,k+Q 2 = e - S - 
£2 cos 20, with <5 = 2/i, 5 2 = eom(m x — m y ) / (2m x m y ) , 
and 9 = tan -1 k y /k x [42]. 

Electrons with spin a of the band i are created by 
the operators cj ka , yielding the non- interacting Hamil- 
tonian: 



Ho = J~] £i,kc| -ka Ci,kq 
i,k 



(1) 



Here and for the rest of the paper the summation over 
repeated spin indices is assumed, and we shift the mo- 
menta of the fermions near the X and Y Fermi pock- 
ets by Qi and Q 2 , respectively, i.e. £x,k+Ch -> £x,k, 

£y,k+Q 2 £y,k- 

As discussed in Ref. [43], this model has eight 
fermionic interactions u n which can be decomposed into 
the spin density- wave (SDW), the charge density- wave 
(CDW) and the pairing channels. Since our goal is 
to study stripe magnetism and the accompanied Ising- 
nematic order, we keep only the interactions in the spin 
channel with momenta near Qi and Q 2 , restricting the 
interacting Hamiltonian to 



1 

"2 Uf 



pin / °2,q 



(2) 



where s iiq = YJ fc cf, k+qQ <x Q( gc; ik( g is the electronic spin 
operator, with Pauli matrices o a p. The coupling u sp ; n 
is the combination of density-density and pair-hopping 
interactions between hole and electron states {u\ and 113 
terms in the notation of Ref. [44]): 



ui t 



-(•••) 



U3 t 



+ (•••) 



(3) 



where the dots stand for the terms with S at pS lt s, which 
only contribute to the CDW channel. Combining the 
two contributions for the SDW channel, we find u sp i n = 
Mi + U3. Once u sp in exceeds some critical value (which 
gets smaller when So and 82 decrease), static magnetic 
susceptibility diverges at (0,7r) and (tt, 0), and the sys- 
tem develops long-range magnetic order. An excitonic- 
type SDW instability in Fe-pnictides, resulting from the 
interaction between hole and electron pockets, has been 
considered by several authors [41, 45-53]. 

Our calculations are done in two steps. In the 
first step, we introduce the bosonic fields A( X ,y) oc 

Sk c r ka' T a^ c (-f ,Y)tep f° r the collective magnetic degrees 
of freedom, integrate out the fermions, and obtain a 
Ginzburg-Landau (GL) action for Ax and Ay. We 
show, in agreement with earlier results [41], that in a 
mean-field approximation only one of the magnetic order 
parameters - either (Ax) or (Ay) - becomes non-zero in 
the magnetically ordered state. This leads to stripe- type 
SDW order in which spins are ordered ferromagnetically 
in one direction and antiferromagnetically in the other, 
i.e. the ordering momentum is either (n, 0) or (0, 7r). 
In the second step, we include fluctuations of the Ax,y 
fields, introduce the collective Ising-nematic bosonic vari- 



able 



A y together with ip cx A 2 X + A r , in- 



tegrate over Ax and Ay, and obtain an effective ac- 
tion in terms of 4> and ip. We analyze this action and 
check whether the system develops an instability towards 
(<j>) ^ before (Ax) or (Ay) becomes non-zero. 



A. The action in terms of Ax and Ay 

A straightforward way to obtain the action in terms 
of (Ax) and (Ay) is to start with the fermionic Hamil- 
tonian W = Wo + "Hint in Eqs. (1) and (2), write the 
partition function as the integral over Grassmann vari- 
ables: 



Z cx 



(4) 



and then decouple the quartic term in fermionic operators 
using the Hubbard-Stratonovich transformation: 



e 2 



where, in our case, 



7+yx) 



Sj.O 



E 



k ^T.ko 



(5) 



r a/3C(x,y),k/3 



and 



y 



HX,Y)- 



We then integrate Eq. (4) over 



fermionic variables using the fact that after the Hubbard- 
Stratonovich transformation the effective action becomes 
quadratic with respect to the fermionic operators. The 
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result of the integration is recast back into the exponent 
and the partition function is expressed as: 



Zoc j dA x dA Y c- s ^ x -^ Y] 



(6) 



If the relevant Ax and Ay are small, which we assume 
to hold even if the magnetic transition is first-order (we 
present the conditions on the parameters below) , one can 
expand S c g [Ax, Ay] in powers of Ax and Ay and ob- 
tain the Ginzburg-Landau type of action for the order 
parameters Ax, Ay. For uniform Aj, the most generic 
form of S cS [Ax, Ay] is 

S cS [Ax, Ay] = r (A 2 X + ± 2 y) + 1 (Ax + A|-) 2 



(A 



x 



A y ; 



v (Ax • A Y )l7) 



Carrying out this procedure, we obtain the coefficients 
r , u, g, and v in terms of the non-interacting fermionic 
propagators convoluted with Pauli matrices (details can 
be found in Appendix A). The coefficient v vanishes in 
our model because of the anti-commutation property of 
the Pauli matrices: o % <ji + crier' 1 = for i ^ j. To 
get a non-zero v, one needs to include direct interac- 
tions between the two electron pockets [41]. The other 
three prefactors are expressed via fermionic propagators 



ro 



u = 



g = 



2 r 

h 2 / Gr,kGx,k 

^spin Jk 

2 J G r,k (Gx,k + G Y .k) 2 
- 2 J G v,k (Gx,k - G Y ,kY 



(8) 



where J k =TJ2 n f (fJ? and k = ( k > w n)> witn momen- 
tum k and Matsubara frequency uj n = (2n + 1) nT. Sim- 
ilar coefficients were found in Ref. [54] , which focused on 
the magnetic instabilities in a two-band model. 

Evaluating the momentum integrals and summing over 
Matsubara frequencies we obtain that J, Gr,kGx,k is 
negative and at perfect nesting diverges as Np log A/T, 
where Np ~ m/(2n) is the density of states of the hole 
pocket and A is the upper cutoff for the low-energy the- 
ory - i.e. the scale at which corrections to the parabolic 
dispersion become of the order one. Away from perfect 
nesting, the logarithm saturates at T = at some fi- 
nite but still large value [55]. As a result, r decreases 
with decreasing T and, if the amplitude of the interaction 
u S pin is not too small, changes sign at some temperature 
T = Tjv,o- This is the mean-field transition temperature 
below which magnetic order appears. Near Tx,o one can 
expand r as r = a(T — Tjv,o), with a > 0. 

The GL action can be straightforwardly extended to 
include the momentum and frequency dependence of the 
A fields, Aj — > Ai(q). This extension does not modify in 



any relevant way the pre-factors of the quartic terms, but 
changes the pre-factor for the quadratic term at small q 
to: 



r a -> Xi,q 



-t-spin 



' / Gr,kGi t k+ q 

Jk 



= r Q +7|z/„| + / q 



(9) 



Here i = X, Y, v n = 2nTn is the bosonic Matsubara 
frequency, 7 is the Landau damping coefficient, and / q is 
in general an anisotropic function of q, / q = q\(\ ± r\) + 
<7y ( 1 =F 77) + r/^qf , with — 1 < r\ < 1 and upper (lower) signs 
referring to Ax (Ay). We verified that all results which 
we obtain below do not depend on whether 77 is finite or 
zero (see Appendix B). To simplify the presentation, we 
then set r\ — and approximate xTq by r + 7|t / „| + <j| + 

i] z ql, the same for Ax and Ay. Below we consider this 
anisotropic 3D dispersion and also an isotropic dispersion 
in dimensions 2 < d < 4, in which case gjj + f] z q\ is 

replaced by the ci-dimensional momentum amplitude q 2 . 

For the two fourth-order terms in (8) we obtain at a fi- 
nite T, expanding around perfect nesting to leading order 
in the chemical potential \i and in the mass anisotropy 
8m = m(m x — m y ) / (2m x m y ): 



7C (3) N F 

,, 0.024. (^) 2 

for Sm < T/e < 1. At T = 0, we obtain 

Np 



(10) 



9 



( e 5m\ 



(11) 



for Sm, <C \fj,\ /so <C 1. In general, u > and u > g. i.e. 
the action (7) increases when either Ax or Ay becomes 
large. If this condition was not satisfied, the expansion of 
S e g in powers of Ax and Ay would only make sense if 
higher order terms were considered. The crucial result for 
our consideration is that g vanishes for circular electron 
Fermi surfaces (Sm = 0) , but is positive for any non-zero 
ellipticity, independent on the sign of eo and on whether 
m x is larger or smaller than m y . 

The action S e s is exact and includes all fluctuations of 
the two bosonic fields. Before we consider these fluctua- 
tions, let us analyze Eq. (7) in the mean-field approxima- 
tion. The conventional way to justify mean-field theory 
is to extend the original one-flavor fermionic model to 
L flavors, such that Ax and Ay couple to L clones of 
Sj, q =o, and take the limit L — >• 00. To do this properly, 



one has to rescale simultaneously u s 



ijL. Af- 



ter this rescaling, the effective action can be written as 
S c s — LS e s, where S e s is the same as Eq. (7), but with 
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rescaled u S pin- Because of the overall factor L the action 
5cff can be approximated by its value at the minimum, 
as corrections to the partition function from fluctuations 
of Aj are small in (logL) /L. 

After solving for the minimum of S c s [Ax, Ay] in Eq. 
(7), we find that, when g = 0, the ground state has a huge 
degeneracy because any configuration A = (Ax) e lC * 1 ' r + 
(Ay)e ,C!l ' r with (Ax) + (Ay) = —r /u minimizes 
S c ff. A non-zero g gives rise to the additional coupling 
2gA x Ay, which breaks this degeneracy. For a positive 
g, this term favors the states in which only one order 
parameter has a nonzero value, i.e. configurations with 
either (Ax) ^ or (Ax) ^ 0, but not both. 

To relate (Ax) or (Ay) to the magnetic ordering, 
we return to the effective action written in terms of 
the double functional integral over the fermionic and 
A fields. Assuming that the electronic spin s ijq=0 = 
c r ka^apCi^p acquires a non-zero magnitude (si,o) 7^ 
0, we again minimize S e s, but this time with (sj,o) as a 
parameter, yielding (s^o) — (A») /w S pin- Since anon-zero 
sx.o (sy,o) implies magnetic order with the momentum 
Qi = (tt, 0) (Q2 = (0, 7t)), the fact that only one of the 
Aj orders means that the magnetic ordering has a par- 
ticular momentum. One can easily verify that in such 
a state spins order ferromagnetically along one direction 
and antiferromagnetically along the other one. 



B. The action in terms of the Ising-nematic order 
parameter 

Since the action (7) is invariant with respect to the 
interchange between Ax and Ay, the onset of either 
(71", 0) or (0, 7r) SDW state breaks not only the conven- 
tional 0(3) spin-rotational symmetry, but also the ad- 
ditional Z2 (Ising) symmetry associated with choosing 
between A x and Ay (Refs. [28, 32, 33]). The issue 
we now consider is whether the Z2 symmetry breaking 
preempts the 0(3) symmetry breaking, i.e. it happens 
before the Ginzburg-Landau parameter ro changes sign 
and the magnetic order sets in. 

Such a Z 2 symmetry breaking without magnetic or- 
der would imply that fluctuations associated with one 
of the bosonic fields are larger than the fluctuations as- 
sociated with the other one, e.g., (A x ) > (A y ) while 
(Ax) = (Ay) — 0. A direct experimental detection of 
this state could be done by performing inelastic neutron 
scattering in detwinned samples and measuring the spec- 
trum at (7T,0) and (0,7t). 

That the action (7) can potentially lead to a preemp- 
tive Ising-nematic instability is evident from the presence 
of the term g (A x — A y ) , which can give rise to an or- 
dered state with (A x ) — (A y ) ^ in a way similar 
to how the Si iq Sj ; _ q term in the Hamiltonian (2) gives 
rise to a state with non-zero (sj,o) 7^ 0. The preemptive 
Ising-nematic instability, however, does not appear in the 
mean-field limit of L — > 00 fermionic flavors simply be- 



cause magnetic fluctuations are absent at L = 00, and a 
non-zero (A?) 7^ appears simultaneously to (Aj) 7^ 0, 
once ro changes sign. However, it may well happen once 
we return to the original model with L = 1 fermionic 
flavor and include magnetic fluctuations. 

To study a potential preemptive Z2 transition, we need 
to introduce collective variables of the fields Ax and Ay. 
Let us introduce auxiliary scalar fields 4> for A x — A y 
and ip for A x + A y . The field tp always has a non- 
zero expectation value (ip) 7^ 0, which describes Gaussian 
corrections to the magnetic susceptibility xiq m Eq. 9. 
Meanwhile, the field (f> may or may not have a non-zero 
expectation value. If it does, it generates a non-zero value 
of (A x — Ay) and the system develops an Ising-nematic 
order. 

The effective action in terms of <j> and ip is obtained by 
using again the Hubbard-Stratonovich transformation of 
Eq. (4), but this time the variable x is either A x + A y or 
A x — A y . After applying this transformation, we express 
the partition function in terms of double integrals over 
the fields (Ax, Ay) and (cf>, ip): 

Z<x J dAxdAy# ( 2V'e~ Soff[AlA ' / ' 1 (12) 

where: 

/*(A 2 X +A|.)+ / *(A 2 Y -Af-) (13) 

Once (p becomes non-zero, we have from (13) 

(Ax) - (Ay) = t (14) 

If the magnetic long-range order is not developed, i.e. 
(Ai) = 0, then it is straightforward to integrate over the 
fields (Ax, Ay). Carrying out the integration we obtain 
the effective action in terms on </> and ip in the form 

s -«[«i=/{C-£ + fM(v + *) 2 -* 2 ]} 

(15) 

We later modify the derivation of S e g [</>, tp] to the case 
where the system has magnetic order. In the next two 
sections we will study the effective action (15) in its 
own mean-field theory and then using a renormalization- 
group (RG) formalism. 

III. MEAN-FIELD THEORY 

To justify a mean-field treatment of S e ff from (15) we 
do a trick similar to the one before, but with the bosonic 
rather than the fermionic variables. Namely, we extend 
the number of components of the fields Ax.y from the 
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original value N = 3 to arbitrary N, assume that the 
<j) and tp fields interact equally with all components of 
Ax.y, and rescale the coupling constants to g — > g/N 
and u -4- u/N. The effective action (15) then has an 
overall pre-factor N implying that for large TV it can be 
analyzed by just taking its value at the extremum. 

We emphasize that the extension to TV — > oo compo- 
nents is a different approximation than the previous one, 
in which we made L — > oo copies of the fermionic fields. 
For L — > oo, the terms of the effective action contain- 
ing the field <j> have an overall pre-factor of 1/L — > 0. 
This means that fluctuations of <f> are large, preventing 
a preemptive Ising-nematic order. On the other hand, 
by abandoning the large L limit and taking instead the 
N — » oo limit, we allow strong fluctuations of both A x 
and Ay , but eliminate fluctuations of the scalar fields <fi 
and tp which only account for a small correction, of or- 
der (logiV) /TV, to the partition function. Rewriting the 
action of Eq. (15) for N components of the Aj field gives: 



S eS [<j>,iP]=N 



2 tp 2 1 

— log 

2g 2u 2 & 



(16) 

The mean-field theory for this action is the saddle- 
point solution of (16), i.e. the minimum of S e ff[<j), tp] 
with respect to <p, and the maximum with respect to tp. 
That e~ S " !! ^'^ increases when tp gets larger is related 
to the fact that Gaussian corrections to ro coming from 
the quartic term (A x + Ay) 2 are confined to the upper 
cutoff of the theory. These corrections just renormalize 
r by a constant, independent on (p 7 and play no role in 
our analysis. What matters to us is how tp is affected 
by <p and vise versa. These mutually influencing terms 
are independent on the upper cutoff and are therefore 
well described within the mean-field theory. Differen- 
tiating Eq. (16) with respect to <p> and tp and taking 
dS c e[(p,tp] /dtp — dS e s [cp, tp] /dtp = 0, we obtain two 
non-linear coupled equations for <p and tp: 



t 
u 



rp + tp + q 2 + tKJ 
q (r + tp + q 2 +7KI) 2 - 



, (ro + ip + q 2 



(17) 



It is convenient to re-express the first formula as a 
self-consistent equation for the renormalized mass of the 
bosonic field, r = r + tp. In the paramagnetic phase, 
r oc £~ 2 , where £ is the magnetic correlation length. It is 
also useful to remove the high-energy contribution to tp 
(i.e. the contribution coming from the upper cutoff), by 
incorporating it into the renormalization of ro. Specifi- 
cally, we rewrite the first equation as: 



r = r + tp = r Q + u 



r + q 2 +7KI 
(r + q 2 +l\v n \f ~ 



(18) 



r + u 



r + q 2 + 7KI _ 1 
(r + g 2 + 7 |^|) 2 -0 2 q 2 +lWn\ 



where ro = ro + u q i_^ \ v 1 ■ For classical systems, only 
the v n = term matters. The remaining momentum 
integral in the last line in (18) is infrared and ultravio- 
let convergent for 2 < d < 4, such that the upper limit 
of integration can be safely extended to infinity. Since 
ro and fo differ only by a constant, fo is also a mono- 
tonic increasing function of T and can be expressed as 
fo = a (T — T/v,o) • For quantum systems such a renor- 
malization of r is not enough in d > 2, and the d + 1 
dimensional integral over momenta and frequency in (18) 
is still generally confined to the upper cutoff. In this sit- 
uation, we use additional renormalizations (see below) to 
restrict the consideration to small energies, at which the 
effective action describes universal low energy behavior. 

For = 0, magnetic order emerges when r = 0, i.e. 
when the static susceptibility for the fields Ax and Ay, 
X^ 1 = r, diverges. The relationship between r and fo 
follows from (18): 



r = r a +u 



I 



1 



1 



r + q 2 +"f\v„\ q 2 +j\v n \ 



(19) 



For classical systems the integrals are infrared diver- 
gent for d < 2, meaning that r never reaches zero - this is 
nothing but the Mermin- Wagner theorem. For d > 2 we 
immediately find from (19) that r and fo vanish simulta- 
neously, i.e. in the absence of a preemptive Ising insta- 
bility long-range magnetic order appears at r = fo = 0. 

Suppose now that a non-zero solution for (f> appears al- 
ready at r > 0, i.e. when the system is still in the para- 
magnetic phase. Once (f> becomes non-zero, the static 
magnetic susceptibilities for the fields Ax and Ay be- 
come non-equivalent: 



XX (q 
XY (q 



0) 
0) 



1 



r — . 
1 



(20) 



Now the magnetic transition occurs when r = \<j>\, i.e. 
at a temperature larger than without Ising-nematic or- 
der. In other words, a preemptive Ising-nematic order 
increases the magnetic transition temperature. 

In what follows, we analyze the phase diagrams result- 
ing from the set of non- linear equations (17) in three dif- 
ferent regimes: the classical regime, where thermal fluc- 
tuations dominate, T7 ^> f (Section III- A); the quan- 
tum regime, T = (Section III-B); and the regime of 
intermediate temperature T7 ~ f , where thermal and 
quantum fluctuations are equally important (Section III- 
C). In all cases, the key parameter that controls the char- 
acters of the magnetic and nematic transitions is the di- 
mensionless ratio: 



u 
9 



(21) 



which measures the strength of the nematic coupling g in 
units of the magnetic coupling constant u. As it is evident 
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from Eqs. (10) and (11), a depends on the parameters 
describing the band structure dispersions. It is therefore 
affected by changes in the chemical potential and in the 
ellipticity - see Section VI for a systematic analysis of a 
as function of doping, pressure, and lattice stiffness. 

In the classical regime (Section III-A) and in the in- 
termediate regime (Section III-C), the other independent 
variable in the phase diagrams (besides a) is the temper- 
ature. In the former case, the temperature dependence 
appears only via the difference from the mean-field Neel 
transition temperature, T — T/v,o- In the latter, we con- 
sider the explicit dependency on T. On the other hand, 
in the T = quantum regime (Section III-B), the inde- 
pendent variable is the distance to the mean-field critical 
point at ro = 0. 

Our primary interest is to obtain the phase diagram 
in the most general case of a quasi-two-dimensional sys- 
tem with anisotropic magnetic interactions. To better 
understand the results in this case, we first consider a 
few limiting cases: 

• We first analyze the classical d = 2 case in Section 
III-A-1. In this situation, there is no finite tem- 
perature magnetic transition, but there is always 
an intermediary paramagnetic phase with nematic 
order. At small a (large nematic coupling g), the 
nematic transition is first-order, whereas at larger 
a (smaller g), it is second-order. The d — 2 case is 
investigated also in the quantum limit (Section III- 
B-l) and in the intermediate-temperature regime 
(Section III-C-1). The behavior in the intermedi- 
ate regime is the same as in the classical regime. 
At T = 0, we find instead that the magnetic and 
the nematic transitions are simultaneous and first- 
order. 

• Next we analyze the classical d = 3 case with 
isotropic magnetic interactions in Section III-A-2. 
We find that, regardless of the value of the dimen- 
sionless coupling a, the magnetic and nematic tran- 
sitions are always simultaneous and first-order. We 
obtain a similar result in the quantum case (Sec- 
tion III-B-2) and in the intermediary temperature 
regime (Section III-C-2). 

• To model the more general case of anisotropic 
quasi-2D systems, we first consider in Section III- 
A-3 the classical system with isotropic dispersion 
in arbitrary dimensionality 2 < d < 3. We find 
that, for large values of a (small g), the system be- 
havior is the same as in the d — 2 case, with split 
second-order magnetic and nematic transitions. On 
the other hand, for small a (large g), the behav- 
ior is similar to the d — 3 case, with simultaneous 
first-order transitions. The two regimes are sep- 
arated by two tri-critical points and an interme- 
diary regime with a first-order nematic transition 
split from a lower-temperature second-order mag- 
netic transition. The quantum and intermediary- 



temperature regimes are analyzed in Sections III- 
B-3 and III-C-3, respectively. 

• In Section III-A-4, we obtain the classical phase 
diagram of a quasi-two-dimensional system with 
anisotropic out-of-plane magnetic dispersion, char- 
acterized by the parameter rj z . For small rj z (weak 
anisotropy) , the phase diagram is the same as in the 
previous case of intermediary dimension 2 < d < 3 
and isotropic interactions. For larger r\ z (moderate 
anisotropy) , we still obtain a regime of split second- 
order nematic and magnetic transitions at large a 
and a regime of simultaneous first-order transitions 
at small a. However, the two regimes are now sep- 
arated by an intermediary regime with a higher- 
temperature second-order nematic transition split 
from a first-order magnetic transition. In Sections 
III-B-4 and III-C-4, we present the results for the 
system behavior in the quantum and intermediary- 
temperature regimes, respectively. 

• Finally, in Section III-A-5, we analyze in more 
details the intermediary regime of split second- 
order nematic and first-order magnetic transitions, 
which appear in the classical phase diagram of the 
quasi-two-dimensional system with moderate out- 
of-plane anisotropy. In particular, we show that the 
first-order magnetic transition is simultaneous to a 
meta-nematic transition, where the nematic order 
parameter jumps between two finite values. 



A. Classical phase diagram as a function of u/g 

When the temperature is high enough (T7 ^> f ), the 
dominant contribution to the sum over Matsubara fre- 
quencies, J q — J d d q/(2w) d in Eqs. (17), comes 
from the term with zero Matsubara frequency. This ap- 
proximation substantially simplifies Eqs. (17) as the re- 
maining momentum integrals can be evaluated exactly. 
We begin with the d = 2 case. 



1. The case d — 2 

For d = 2, the integration over momentum in Eqs. (17) 
yields: 

r = fo - - log (r 2 - cj) 2 ) 

r = ^coth^y^j (22) 

where we defined the renormalized parameters f — 
r + SlogA, u = uT/(2n), and g = gT/(2n), with A 
denoting the upper momentum cutoff. The second equa- 
tion implies that a solution with <p 7^ is only possible 
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Figure 3: (Color online) fo as the function of the Ising- 
nematic order parameter 0* = 2<j>/g for two representative 
values 1 < a < 2 (blue/dark line) and a > 2 (orange/light 
line). The value of (f>* at which the solution first appears 
upon decreasing fo sets the type of the nematic transition. 
For a > 2, the solution first emerges at 0* = 0, and the tran- 
sition is second-order. For 1 < a < 2, fo is non-monotonic 
function of cj>* and the solution first emerges at a non-zero 0* . 
The nematic transition then becomes first-order. 



when r > 0. Eliminating r from these equations, we 
obtain after further rescaling: 



coth < 



a log 



( <t>* 



V sinh ( 



(23) 



2r /g - a log (5/2), and a = u/g 



where (/)* = 2<f>/g, ro 
(see Eq. 21). 

Recall that the original model is constrained to u > g, 
i.e. a > 1. The variable f decreases with decreasing 
T, since it only differs by a constant and by an overall 
factor from ro. The leading instability of the system upon 
decreasing T is into a state with 4>* corresponding to the 
maximum of the left-hand side of Eq. (23). A simple 
analysis shows that the maximum is at <jf — for a > 2 
and at a finite 4>* for 1 < a < 2 (see Fig. 3). The 
implication is that, for a > 2, the system undergoes a 
second-order Ising-nematic transition at fg r = 1 (r = 
g/2), while for 1 < a < 2 the Ising-nematic transition 
becomes first-order and the solution for <jf first appears 
at a larger f = f™ ax (see Figs. 3 and 4). The value 4>* r 
at which the left-hand side of Eq. (23) has a maximum 
gradually increases as a decreases, approaching infinity 
as a — > 1. 

The actual f^ r at which the first-order Ising-nematic 
transition occurs is somewhat smaller than f™ ax because 
at f™ ax the effective action S e g only develops a local 
minimum at a non-zero <j)* r . The actual transition oc- 
curs when the value of the effective action at the local 
minimum becomes equal to S c g((j)* = 0). To obtain fp r 
we then need to evaluate the effective action at both min- 



ima 0* 7^ and 



0, and find fff at which the two 



terms become equal. For better clarity, we compute S e s 
not only at the extrema (when the self-consistent equa- 
tions (22) are valid), but for arbitrary <jf at a given fo- 



To do this, we solve the equation dS c g/dip = to obtain 
r (a, </)*,fo), substitute it back into the effective action 
and obtain S c s [a, <f>* , ro]- Carrying out the calculations, 
we obtain, neglecting a constant term, 



Sett [a, 



,ro. 



9 S 
— - Jeff 



a, 



,r \ 



(24) 



where 



Sett 



-log 

4 S 



and r* = r* 



(r*) 2 - on 2 



r* 2 -log 



log 



(25) 



, a, ro) is the solution of the equation 



a 

r + — log 



(r*) 2 -W>T 



(26) 



In Fig. 4 we plot S e g[a, 0*,r o ] for a = 1.5 as a 
function of <f>* for several values of fo. We see that 
for fo > 1.21, S e g has a minimum only at <jf = 0, 
and Seg [1-5, </>* , fo] monotonically increases with increas- 



ing 



However, once ro becomes smaller than r Q 



which for this value of a is f™ ax = 1.21, the function 
Sea [1.5, (j>* , fo] develops an inflection point at <jf cr f=a 2.72. 
At smaller fo, this inflection point gradually splits into 
a maximum at d> < d>* and a minimum at d> > d>* . 



At 



= 1.15 the values of S e s at </>* = and at the 



local minimum <\>* ^ become equal, and the system un- 
dergoes a first-order Ising-nematic transition. The local 



minimum of S c s at <j)* = survives down to 



= 1. 



Below this temperature, the effective action has only one 
minimum at a finite </>*. 

The jump in <f> at the first-order transition affects the 
susceptibilities associated with the two magnetic order 
parameters and Ay, which become non-equivalent 
once <p becomes finite, see Eq. (20). This implies that the 
static susceptibility and the magnetic correlation length 
change discontinuously at the first-order Ising transition, 
even though there is no magnetic instability (the value to 
which <f) jumps is always smaller than r, see the second 
equation in (22)). Actually, as we already mentioned, 
magnetic order never emerges for d = 2 at a finite tem- 
perature, so Z2 is the only symmetry that gets broken. 

To obtain the phase diagram in the (a, fo) plane, we 
need to analyze the behavior of the system once the Ising 
order sets in. We found that the Ising-nematic order pa- 
rameter continuously increases with decreasing fo for all 
a > 1, implying that there is no other first-order transi- 
tion line in the phase diagram besides the one at which 
the Ising-nematic order develops. The phase diagram is 
presented in Fig. 5. The upper spinodal in this figure 



corresponds to 



where a local minimum of S e s ap- 



pears at <p* 7^ 0, and the lower spinodal refers to f™™ 
where <jf — ceases to be a local minimum of S e e- A 
first-order transition happens in between the upper and 
lower spinodal lines. 
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Figure 4: (Color online) The effective action S e a [a, 0*,fo] 
from Eq. 25 as a function of <j>* for a — 1.5 and several fo. 
The evolution of S c a is typical of a first-order phase transition. 
At large fo, S e ff [1-5, </>*,fo] monotonically increases with in- 
creasing 0*. At f™ ax = 1.21, the effective action develops an 
inflection point (black/dark solid curve) which at smaller fo 
splits into a maximum and a minimum. At fg r = 1.15, the 
values of S c s at <jf>* = and at the local minimum at finite <j)* 
become equal (black dotted curve), and the system undergoes 
a first-order Ising-nematic transition. The local minimum in 
S e ff at (f>* — survives down to f™ ln — 1 (blue/light gray 
dotted curve). 
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Figure 5: (Color online) Calculated phase diagram in the 
(a, fo) plane for the d = 2 case and high enough temperatures, 
where classical fluctuations dominate. In the region a < 2, 
the nematic transition is first-order. The dashed blue/dark 
line represents the upper spinodal (where the ordered state 
becomes metastable), while the dotted blue/dark line denotes 
the lower spinodal (where the disordered state ceases to be 
metastable), and the dashed-dotted orange/light line is the 
actual transition line where the global minimum of S e s shifts 
from <j> = to <f> ^ 0. We recall there is no magnetic order at 
a finite temperature in d = 2. In this figure, we fixed u = 1 
and changed g. 



2. The case d = 3 

A very different phase diagram emerges in three di- 
mensions. The momentum integration in Eqs. (17) now 
yields 



(27) 



where fo = ro — 2A-u/7r, and, as before, u = uT j (27r), and 
g = gT ' j (27r). It is convenient to re-express these equa- 
tions in terms of r and z = (f>/r (0 < z < 1). Eliminating 
r from (27) we then obtain: 



1 



i + VT^z 1 



(28) 



The right-hand side of (28) is an increasing function 
of z for all values of a. Therefore, the first instability 
of the system is at fo = g 2 (1 + a) /8 towards the state 
with the largest possible value of z, namely, z = 1. As a 
result, the order parameter (f> jumps at the Ising-nematic 
transition to <p = r. For such </>, the static susceptibility 
associated with one of the magnetic order parameters di- 
verges, see Eqs. (20), implying that the first-order Ising- 
nematic transition instantaneously brings the system to 
the verge of magnetic order. The divergence of the sus- 
ceptibility reflects the fact that the Goldstone modes of 



the magnetically ordered state are gapless. In the large- 
N limit, the impact of longitudinal fluctuations, where 
the susceptibility remains finite, is negligible. 

In this situation, we have to extend our analysis and 
investigate the possibility that Ax jumps to a non-zero 
value at the Ising transition. To do this, we go back 
to the effective action (12) written in terms of double 
integrals over both Ax,y and their collective variables cf> 
and ip, and allow one of components of Ax to have a 
non-zero mean value M — (A^-). Expressing A x = 

MSqfiS 1 '^ + A x and integrating over Ax and Ay we 
obtain: 



Seff [<j>,1p,M] =S eS 



M 2 (r 



1 + (r - 



1 



q 2 (r — (f> + q 2 ) 



(29) 



In the spirit of the N — > oo approximation, we rescale 
M 2 -> M 2 N and differentiate S c s [0, ip, M] over all three 
variables. This yields the set of three self-consistent equa- 
tions for </>, r, and M. The equation for M is: 



M(r 



l + (r- 



1 



q q 2 (r-cp + q 2 ) 



0. (30) 



The solution is either M — or r — <f>. We take M ^ 
and r — cj). The other two equations for r — <j> are 



'0 
9 



2(f) + gM l 



(31) 
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where M - 
we obtain: 



MyVm/T. Solving the first equation for 




(32) 



Substituting this solution into the second equation 
yields f as a function of M: 




(33) 



Following the same strategy as before, we look for the 
value of M correspondent to the largest fo, i.e. we de- 
termine the value of M that emerges at the largest tem- 
perature. A straightforward calculations shows that it is 
non-zero and equals to 



M 



4(a-l)' 



(34) 



This implies that the first-order Ising-nematic transi- 
tion at which <fi reaches its maximum value (= r) triggers 
a first-order magnetic transition into a state with a finite 
M = (Ax)- We emphas ize, however, that in this and 
similar cases which we consider below, the magnitude 
of the nematic order parameter <f> is larger than gM , 
i.e. the magnetic transition is secondary to the nematic 
transition. Indeed, substituting M from (34) into (32) 
we obtain 



= a (gM 2 ) > gM 2 



(35) 



We emphasize that, without nematic instability, a 
"pure" magnetic transition would take place at a smaller 
temperature, when the condition r — is satisfied, in- 
stead of r = <j). Therefore, even though both transi- 
tions are simultaneous and first-order, the nematic one 
is the primary transition, and the magnetic transition is 
induced by the feedback from the nematic order. 

Note also that in the mean-field theory the suscep- 
tibility associated with the field Ax remains massless 



(X 



A' 



■ — <fi = 0) despite that magnetic order develops 
and enforces a gap in the spectrum of longitudinal fluctu- 
ations. This is a consequence of taking the N — > oo limit, 
in which one longitudinal mode is negligible compared to 
N — 1 gapless transverse modes. 

Like in the d = 2 case, the value of fo, at which the 
effective action develops an inflection point and the so- 
lution with <j) = r and M ^ first appears, is the upper 
spinodal f ™ ax . The actual first-order transition occurs at 
a smaller f^ r , at which the values of S e s at <j> = M = 
and at the local minimum at (f> = r, M ^ become equal 



d=3 



a 



Figure 6: (Color online) Classical phase diagram in the (a, fo) 
plane for the d = 3 case. For any value of a, the nematic and 
magnetic transition are simultaneous and first-order. 



(for r — <j> the local minimum is with respect to variations 
of M). 

The phase diagram for d = 3 is presented in Fig. 6. 
As in the d — 2 case, there are no additional transition 
lines, i.e. after jumping to finite values at the first-order 
transition, both cj> and M continuously increase with de- 
creasing fo ■ The similarities between the d = 2 and d = 3 
cases, however, stop here because the phases below the 
critical line are very different. In d = 3, Ising-nematic 
order immediately triggers a magnetic order, such that 
there is no regime where the Z2 symmetry is broken but 
the 0(3) symmetry is unbroken. In d = 2, there is no 
magnetic order below the Ising-nematic transition, which 
can be either first-order or second-order, depending on a. 

We now consider what happens in anisotropic three 
dimensional systems. There are several ways to model 
the anisotropy: one can either consider the dimension d 
to be an arbitrary number between 2 and 3 or one can 
keep d = 3 but consider an anisotropic magnetic disper- 
sion with different stiffness along q z and along q Xl q y . It 
turns out that the system behavior is universal at small 
and large a, but at intermediate a it depends on the 
choice of the model. This will lead us to the two phase 
diagrams shown in Fig. 1. We first consider arbitrary 
2 < d < 3 and then an anisotropic d — 3 dispersion. 



3. The case 2 < d < 3. 

We assume first that M = and later extend the 
formalism to include a non-zero magnetization. For 
2 < d < 3, the integration over momentum in Eqs. (17) 
yields 



r °~4 



(r + . 



(r + , 



2 + (r- 



(36) 
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where: 



u 
u 



9 



2TS d 
(2nf 



~d-3 



dx ■ 



1 

T 



(37) 



2 d -%^r 



(f ) 



and S d = 27r d / 2 /r(d/2) is the area of a d-dimensional 
sphere with unit radius. Introducing, as before, z = <f>/r 
(Q < z < 1), we solve the second equation for r (z), sub- 
stitute the solution into the first equation, and obtain: 



with: 



'(*,<*) 



r ° = (f ) ^ ' 



(l + z)~ - (1-z)" 



(38) 



d-2 
4-d 



d-2 / 1 \ d-2 / | 

(l + z) 2 ( a+ - j + (1 - z) 2 (a-~ 



(39) 



We determine the leading instability of the system by 
again looking at the location of the maximum of Q (z, a) 
at a given a. In Fig. 7, we plot Q(z,a) for d = 2.5 and 
three different values of a. We find three different regimes 
for the behavior of the system: For 1 < a < a c i = 2, the 
first instability is at z = 1; for a c \ < a < a C 2 = 3.5, the 
first instability is at an intermediate value < z < 1; 
for a > a C 2, the first instability is at z — 0. Expanding 
Q (z, a) near z — and z = 1 and analyzing the sign of 
the slope, we find the expressions for a c \ and a C 2 in an 
arbitrary dimension: 



3-d 



a C 2 



6-d 
6 -2d' 



(40) 



For d — 2 we have a c \ = 1 and a C 2 = 2, i.e. the 
regime 1 < a < a c \ disappears, in agreement with Fig. 
5. On the other hand, for d = 3, both u\ and «2 tend to 
infinity, and the region 1 < a < a c \ extends to all values 
of a, in agreement with Fig. 6. For d between 2 and 3 all 
three regions are present, as shown in the phase diagram 
of Fig. 8. 

In the first region 1 < a < a c \ (region I in Fig. 8), the 
behavior of the system is the same as we found in d = 3: 
(f> jumps at the Ising-nematic transition to the largest 
possible value <f> = r, triggering a simultaneous magnetic 
transition. To determine whether the latter is also first 
order, we extend the analysis of the effective action in 
the same way as we did for d = 3: we introduce magnetic 
long-range order via the order parameter M = (^x) an< ^ 
obtain the set of three coupled equations for M , <fi, and ip. 
The equation for M again gives either M — or <p = r. 
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Figure 7: (Color online) Q (z, a) oc fo, as denned in Eqs. 
(38) and (39), plotted as function of z = <j>/r for three rep- 
resentative values of a at d = 2.5: (a) a > a C 2 = 3.5, (b) 
2 = a c i < a < a C 2, and (c) a < a c i- The maximum of Q 
shifts from z = at large a to intermediate z at intermediate 
a, and to z = 1 at small a. 



We choose r = 
as 



and re-express the other two equations 



ro 
fj 



uM l 



(20)- 



gM 2 



(41) 



with M/M_= yjuju (see Eq. 37). For a < a cl , we 
expand in M and obtain an explicit equation relating fo 
to the magnetization M: 



ro 



g\ 4-d / 1 + a 



(I) 



2(4 -d) 



Q(M,a), (42) 
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Figure 8: (Color online) The calculated phase diagram in the 
(a, fo) plane for the representative case d = 2.5. In region I, 
there is a simultaneous nematic and magnetic first-order tran- 
sition. In region II, the two transitions split, but the nematic 
transition remains first-order, while the magnetic transition 
becomes second-order. In region III, the nematic transition 
also becomes second-order. The three regions are separated 
by two tri-critical points. The temperature behavior of the 
nematic and magnetic order parameter in each region is shown 
in Fig. 9 below. 



with: 



f)>/(t> ,M <l>/<t>o,M 




(c) (d) 



Figure 9: (Color online) Nematic order parameter <j> 
(blue/dark gray curves) and magnetic order parameter M 
(red/light gray curves) as function of — Afo = ra tCr — ro oc 
(T s — T) for four different values of a at d = 2.5. In (a) and 
(b) we use a > a C 2 = 3.5 (region III of the phase diagram), 
in (c) we use a C 2 > a > a c \ = 2 (region II of the phase 
diagram), and in (d) we use a < a c \ (region I of the phase 
diagram). (j>o is the value of <j> corresponding to <fi = r. 



Q (M, a) « 4 (3 - d) (a cl - a) M 2 



(I) 



For a < ct ci the maximum of Q (M, aj is at a finite 
magnetization: 



M 



(4-rf) (3-d) (2gy 
(l + a) (d-2) 



1/2 



(44) 



Therefore, in region I, the first-order Ising transition 
triggers the first-order magnetic transition. We again em- 
phasize, however, that the Ising-nematic order parameter 
is larger than the square of the magnetic order parame- 
ter, i.e. the magnetic transition is the secondary transi- 
tion, triggered by the preemptive Ising transition. This 
is most clearly seen from Eq. (44): the jump in the mag- 
netization approaches zero as a approaches a c ±, while 
the jump in the Ising-nematic order parameter remains 
finite and reaches <j) = (l/2)(g/2) 2 /^ 4 ' d \ Thus, a cl cor- 
responds to a magnetic tri-critical point, while a preemp- 
tive first-order Ising transition exists on both sides of it. 
In Fig. 9 d, we present numerical results for cf> and M as 
functions of fo in region I. 

In the second region a c ± < a < a C 2 (region II in Fig. 
8), the Ising-nematic transition is still first-order, but 



the magnitude of the jump of (j) is smaller than the value 
required to trigger a magnetic transition, i.e. <j) < r. As 
a result, in this region the magnetic and Ising transitions 
are split, with the former occurring at a smaller fo (i.e. 
at a smaller temperature) than the latter. The magnetic 
transition then becomes second-order, as the maximum 
of Q (M, a) remains at M = for a > a c i, see Eq. (43). 
In Fig. 9 c, we present the numerical solution for <f> and 
M in this region, as function of temperature. 

To determine what happens as the system approaches 
a C 2, we expand Q(z, a) in Eq. (39) around z = and 
obtain: 



Q{z,a) 
(d-2)^ 



2a 
d-2 



1 



(3-d) 



. , 2 (15 + 2d)(6-d) , , 4 

(cv C 2 — a) z H (Kct C 2 — a) z 



with a constant k 



120 



90-7<i-tr 
(6-d)(15+2d) 



45) 



< 1. Clearly, the 



jump of the Ising-nematic order parameter across the 
transition approaches zero as a increases towards a C 2, 
Zmax ~ V a c2 — ol. Therefore, a C 2 is a nematic tri- 
critical point, beyond which the Ising-nematic transition 
becomes second-order (region III in Fig. 8). In this re- 
gion, a non-zero <j) gradually develops below the transi- 
tion line, and the magnetic transition splits even further 
from the Ising-nematic transition. We present the nu- 
merical solution for and M in region III in Fig. 9 b,a. 



14 



4>/r 



Ar s P lil 



1.0 



0.5 



a = 3.6 



a = 7 



-Ar n 



0.02 



0.04 



0.06 



Figure 10: (Color online) Nematic order parameter <f>, in units 
of the renormalized magnetic mass term r, as function of the 
reduced temperature — Afo = ro, C r — ro oc (T s — T) at d = 2.5 
for two different values of a : a = 7, which is far from the 
nematic tri-critical point a C 2 = 3.5 (light-green/ light gray 
curve), and a = 3.6, which is very close to the nematic tri- 
critical point (dark-green/dark gray curve). In both cases, <f> 
evolves continuously from zero, but the slope is very large for 
a only slightly above a c i- 



Note that the slope with which (j> increases below the 
instability remains very high over some range of a in 
region III, as illustrated in Fig. 9(b) and, more trans- 
parently, in Fig. 10. From the practical point of view, 
the Ising transition can then still be viewed as almost 
first-order. This can also be seen from Eq. 45 because 
in coefficient of the z 4 term, k remains very close to one 
for all 2 < d < 3 (0.048 < 1 - k < 0.054). As a re- 
sult, if a is not too far from a C 2, 4> increases very rapidly. 
For instance, the value fo/rfi 1 ^ for which <j> increases to 
(f) = r/2 is \fo/r^ ax — 1| < 10~ 3 (a — a c2 ), almost inde- 
pendent on the dimensionality. 

We emphasize once again that for a < a C 2, when the 
nematic transition is first order, the value of ro at which 
the solution of Eq. (38) first appears is the upper spin- 
odal of the system, f™ ax . The lower spinodal f™ m is the 
value of fo below which <fi = is no longer the local mini- 
mum of S c fi- The actual phase transition in the presence 
of fluctuations is between f™ ax and f™ m , and is close to 
if fluctuations are weak. If we use f™" 1 instead of 



at which the first-order tran- 



'o 

zvnax 

sition occurs, we find that a C 2 remains intact, but a c \ 
moves to a larger value a' cl given by: 



(46) 



2 

d-2 



One can easily check that for 2 < d < 3, a c \ < 
a' cl < a C 2- Thus, the magnetic tri-critical point does 
not merge with the nematic tri-critical point, even if 
the Ising-nematic transition occurs at the lowest possi- 
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Figure 11: (Color online) Splitting between the nematic and 
magnetic transitions, AfQ Pllt = fQ em — r™ as , as function of 
a = u/g, for the case d = 2.5. The dashed line refers to 
regions I and II of the phase diagram of Fig. 8, whereas the 
solid line refers to region III. 



Notice that in the phase diagram of Fig. 8 the split- 
ting between the magnetic and nematic transitions ini- 
tially increases with increasing a, once the system crosses 
the nematic tri-critical point. However, for small values 
of the nematic coupling g (large values of a) , this split- 
ting is expected to decrease, since it should tend to zero 
for infinitesimal g. This non-monotonic behavior is a 
consequence of the fact that, at smaller a, the strong 
first-order nematic transition brings in a simultaneous 
magnetic transition, despite the fact that the magnetic 
correlation length is rather small immediately above the 
nematic transition. This correlation length, meanwhile, 
increases monotonically as a increases. We show the non- 
monotonic behavior of the splitting between the magnetic 
and nematic transitions as function of a in Fig. 11 for 
the case d — 2.5. Note that the splitting begins to slowly 
decrease only at a > 10. 



4- Anisotropic d — 3 case 

We now show that the phase diagram obtained for the 
case 2 < d < 3 is qualitatively the same as the phase dia- 
gram of the d — 3 model with strongly anisotropic, quasi- 
2D magnetic dispersion. However, for less anisotropic 
systems, we find a different phase diagram, with new be- 
havior at intermediate a. 

To model the anisotropy, we consider the system com- 
posed of stacked two-dimensional layers. In each layer 
n, we define the magnetic order parameters A™, with 
i = X.Y. The coupling between different layers adds the 
term -^A" • A" +1 to the action with r) z < 1. As a 
result, the magnetic susceptibility in Eq. (9) becomes: 



rjz sin 



2 Qz 



(47) 



ble f^™. 



with < q z < 2-7T and qn = (q X T<ly)- An alternative 
possibility is to consider anisotropic but still quadratic 
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dispersion 



X 



1 



ro 



2„2 



(48) 



with < P < 1, and set the same cutoff A for all three 
components of q. We found that the phase diagram is 
the same regardless of whether we use Eq. (47) or (48) 
for the bosonic susceptibility. For definiteness, here we 
present the results for \q given by Eq. (47) and consider 
the case of Eq. (48) in Appendix C. 

We plug Eq. (47) into the self-consistent equations (17) 
and evaluate the three-dimensional integrals. Absorbing 
the cutoff A into ro, fo = ro + 2-0 log (2A), we obtain in 
the paramagnetic phase: 



r - u log (y/r + 4> + \Jr] z + r + qjJ 
-■Slog (y/r - (j>+ 



= 9 log 



y/r + 4 1 + VVz + r + ' 



y/r- <fi + \fr\ z + r- 



(49) 



where, as before, u — uT/ (2ir) and g = gT / (2ir). >From 
the second equation, we obtain r as function of <f>: 



r{cj>) 



' 2 



cosh 



\ 



4 



sinh 2 * 

9 



(50) 



Substituting (50) into the first equation in (49) we ob- 
tain fo (<f>), whose maximum determines the first insta- 
bility of the system. At a non-zero M we find, as before, 
r = (j) and 

ro (<f>) = 

<j) + u [log (v 7 ^! + VVz + 2<fe) - M 2 ] (51) 
where the magnetization is given by: 



M 



fJ 



- log 



2(f> + VVz + 20 



(52) 



Notice the characteristic logarithmic dependence on 
the anisotropy, log?7 z . Using these expressions, we obtain 
the phase diagram of the anisotropic system by analyzing 
the first instability upon decreasing fo. 

We show our results for various r\ z in Fig. 12, where we 
present four possible phase diagrams together with their 
respective fo(z, a) profiles for different values of a (recall 
that z — 4>/r). The phase diagrams shown in panels (a) 
and (d) exist over sizable ranges of small and moderate 
r] z , respectively. On the other hand, the phase diagrams 
in panels (b) and (c) exist only over rather narrow ranges 
of rj z and just show how the system actually evolves from 
the regime (a) into the regime (d). 

The behavior at small rj z shown in Fig. 12 (a) is the same 
as for a fractional dimension 2 < d < 3, Fig. 8. Namely, 
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Figure 12: (Color online) Left panels: the four different phase 
diagrams for the anisotropic 3D model with the magnetic 
susceptibility given by Eq. (47). The behavior found in 
isotropic systems for a fractional 2 < d < 3 only holds at 
small fj z = rj z /g. At larger g z , the magnetic and nematic tri- 
critical points q c i and a c i interchange and a new behavior 
emerges at intermediate a. The behavior presented in panels 
(a) and (d) takes place over sizable ranges of small and mod- 
erate f] z , respectively. The behavior presented in panels (b) 
and (c) takes place only over a narrow range of rj z , around 
the f} z value for which a c i and a c i cross. Right panels: the 
behavior of fo(z) for several different a for each phase dia- 
gram (z = 4>/r). The key difference between the regimes (a) 
and (d) is that in (a) fo(z) has only one maximum at every 
a, while in (d) fo(z) has two distinct maxima and a minimum 
in between. 



there is a magnetic tri-critical point at a c i and a ne- 
matic tri-critical point at a C 2 > a c i, with the regime of 
split first-order nematic and second-order magnetic tran- 
sition in between. The behavior at larger r\ z is, however, 
different. We see from Fig. 12(d) that, in the regime of 
intermediate a, there is a second-order nematic followed 
by a first-order magnetic transition. 

This change in the system behavior upon increasing rj z 
can be better understood by considering how the ratio 
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Figure 13: (Color online) Magnetic and nematic tri-critical 
points « c i (red/light gray dotted line) and a C 2 (blue/dark 
gray solid line) as functions of fj z for the 3D anisotropic model 
with the magnetic susceptibility given by Eq. (47). Notice 
that a c i and a C 2 cross at fj z ~ 0.43. 



of the two tri-critical points a c \ / ol c i evolves with r\ z . In 
2D, when rj z = 0, we have a c ijoi c \ — 2. For the isotropic 
system in 2 < d < 3, we have from (40) that a C 2/a ci = 
3 — d/2, implying that a C 2 remains larger than a c \ for all 
2 < d < 3, what gives rise to the phase diagram of Fig. 8. 
However, there is no requirement that a C 2 must remain 
larger than a c \ as the system approaches the isotropic 
3D regime - in fact, the only requirement is that in this 
limiting case both a c \ and a C 2 diverge. In Fig. 13 we plot 
a c \ and a C 2 as a function of r\ z for the 3D system with 
anisotropic dispersion. We see that a c \ and a C 2 cross at 
a certain ry 2 , beyond which the nematic tri-critical point 
occurs at a smaller a than the magnetic tri-critical point. 
This leads to the phase diagram of Fig. 12(d). 

We can show analytically that a c \ and a C 2 cross upon 
increasing rj z . To obtain a C 2, we substitute r{(f>) from 
(50) in (49) and expand the right-hand side in powers 
of (j). Since the quadratic <fi 2 term vanishes at a C 2, we 
obtain: 



a, = 2+5|l 



(53) 



where fj z = rj z /g. To obtain a c i, we express <f> in terms 
of to using Eq. (52), substitute the result in Eq. (51) 
and expand the right-hand side in powers of M. Since 
the quadratic M 2 term vanishes at a c i, we find: 



y/2^y/fj z + 20. 



(54) 



cl 



where 



is the solution of 



= log 



cl 



'rjz 



(55) 



A simple analysis then shows that a C 2 becomes larger 
than a c i at fj z = 0.43. One has to be careful to prop- 
erly determine a C 2 for fj z > rf z , since the actual nematic 




Figure 14: (Color online) Calculated phase diagram in the 
(a, fo) plane for moderately anisotropic d — 3 system with 
?7 Z = 0.3u. Three different regions of system behavior are 
marked in the same way as in Fig. 8. Notice that the positions 
of the tri-critical points are reversed with respect to the case 
of strongly anisotropic systems. 



tri-critical point takes place when z = is the global 
maximum of tqIz), and not only a local maximum. This 
subtlety does not affect the result that a c \ and a C 2 cross, 
nor the regime fj z < fj z . In the phase diagram shown in 
Fig. 14, a C 2 is the actual nematic tri-critical point. 

The phase diagram shown in Fig. 12(d) contains a new 
type of a phase transition not seen in quasi-two dimen- 
sional systems, namely, a meta-nematic transition. We 
explore it in more detail in the next subsection. 



5. Meta-nematic transition 

In Fig. 14 we explicitly show the three different 
regimes that appear in Fig. 12(d), as we did in Fig. 8. 
The behavior of the nematic and the magnetic order pa- 
rameters as functions of the distance to the transitions in 
all three regimes is shown in Fig. 15. For these particular 
figures, for convenience, we considered rj z — Q.Zu instead 
of r\ z = 0.7<7, which does not change the properties of the 
phase diagram. 

A careful analysis of the profile of tq{z, a) as a function 
of z for the phase diagram in Fig. 12 (d) shows that, 
at small a, the maximum of fo(z,a) is at z = 1 (i.e. 
at cf> = r). This implies that upon decreasing fo the 
system undergoes a first-order nematic transition, which 
triggers a simultaneous first-order magnetic transition. 
This system behavior takes place in region I of Fig. 14, 
and is displayed in Fig. 15 (c). 

However, as a increases, fo(z, a) develops a local max- 
imum at z — as well. At the nematic tri-critical point 
0^2, ro(0, a) becomes equal to fo(l,a), and ?o(z, a) has 
a minimum for < z < 1. Once a becomes larger than 
a C 2, the absolute maximum of fo(z, a) shifts to z = 0, 
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Figure 15: (Color online) Nematic order parameter <f> 
(blue/dark gray curves) and magnetic order parameter m 
(red/light gray curves) as function of — Afo = ro, C r — ro oc 
(T s — T) for three different values of a in the case of a moder- 
ately anisotropic system with r\ z — 0.3m (see Fig. 14). Panel 
(a) is for a > a c \ « 2.91 (region III of the phase diagram), 
panel (b) for a C 2 ~ 2.54 < a < a c i (region II of the phase di- 
agram), and panel (c) is for a < a C 2 (region I of the phase di- 
agram). In panel (b), the dotted line marks the meta-nematic 
transition, coincident with the magnetic transition. <f>o is the 
value of (f> corresponding to <j> = r. 



i.e. the nematic transition becomes second-order. As ro 
decreases, a solution with small and finite z ^ devel- 
ops. But because the profile of r~o(z, a) has two maxima, 
another solution with z = 1 necessarily appears once ro 
gets smaller than some critical value. We checked that 
the value of S c s at z = 1 eventually becomes larger than 
the value corresponding to the solution at small z. As 
a result, at some fh, the nematic order parameter un- 
dergoes a meta-nematic transition, where it jumps from 
some small value ^ to <fi = r [z = 1). In accordance 
to what we found earlier, we explicitly confirmed that 
such a transition triggers a first-order magnetic transi- 
tion. Therefore, for a > a C 2, the second-order nematic 
and meta-nematic transitions split. The second-order 
transition occurs first, and the meta-nematic transition 
occurs at a smaller ro, triggering a first-order magnetic 
transition. This system behavior takes place in region II 
of the phase diagram of Fig. 14, as shown in Fig. 15 (b). 

As a increases further, the magnetic transition remains 
first-order up to a — a c \ > a C 2- At this point, the posi- 
tion of the maximum of ro(M, z = 1, a) shifts to M = 0, 
i.e the magnetic transition becomes second order. We 
found both analytically and numerically that at exactly 
the same a = a c \, the local minimum of fo(z,a) moves 
to z = 1, meaning that fo(z, a) becomes a monotonically 



decreasing function of z for all < z < 1. As a result, 
for a > a c i, the nematic order parameter <j) monoton- 
ically increases with decreasing fo until it reaches the 
value (j) = r, where the system undergoes a second-order 
magnetic transition. This system behavior corresponds 
to region III in Fig. 14, and is the same as region III 
in Fig. 8. Notice from Fig. 15 (a), however, that for 
a gl a c i, the proximity to the meta-nematic transition 
line produces a kink in the temperature dependence of 
the nematic order parameter. 

We acknowledge that, for values of rj z that are not very 
small - such as the one considered here - the magnetic sus- 
ceptibility of the anisotropic 3D system may have extra 
terms not contemplated by our formalism, such as cor- 
rections to the continuous in-plane dispersion. However, 
as we show in details in Appendix C, the same behavior 
is obtained for an anisotropic quadratic dispersion with 
equal momentum cutoff along all three momentum direc- 
tions. This gives extra confidence that the phase diagram 
of Fig. 14 may be realized in at least some moderately 
anisotropic systems. 

We point out that a phase diagram with the interme- 
diate regime of a second-order nematic transition and 
first-order magnetic transition was earlier obtained in a 
semi-phenomenological model for the interaction between 
the structural and magnetic degrees of freedom [56]. Kim 
et al considered a microscopic version of that model [15], 
showing that anharmonic elastic terms may bring the 
system into the regime of split second-order nematic and 
first-order magnetic transitions (we discuss the coupling 
between the lattice and the nematic degrees of freedom 
in Section V below). Our results show that such a behav- 
ior can be obtained in a purely magnetic model, even if 
the coupling to structural degrees of freedom is negligibly 
small. 



B. Quantum phase diagram as a function of u/g 

To complement our analysis of thermal fluctuations, 
we now consider the opposite limit of T = 0, when the 
Matsubara frequency become a continuous variable and 
J = j d d qdv n j '(2n)( d+1 \ Now ro is a function of some 
control parameter, e.g. doping, pressure, or applied field. 
We show that no new phases appear in the T = limit, 
compared to the three phases that we found previously 
in the classical limit. As before, we consider first d = 2, 
then d = 3, and then arbitrary d between 2 and 3 and 
anisotropic 3D systems. 



1. d = 2 

Integrating the self-consistent equations for <f> and r in 
(17) over both v n and q, we obtain for d = 2 

f, A2 r , i <A 

1 = g log - + 1 — - tanh 1 - 

V V r2 - 4> 2 4> r J 
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Figure 16: (Color online) The function fi,{z) from Eq.(58). 



■log 



A 2 



)tanh" 1 ^ (56) 



where g = gj (4ir 2 j^, u = uj (47r 2 7) are dimensionless 
couplings, and A is the upper limit of the integral over 
momentum. The parameter fo again decreases as the sys- 
tem approaches the magnetic transition. Note that the 
dependence on A is still present even after we absorbed 
the r— and ^-independent contribution of the right-hand 
side of Eq. (18) into fo. This dependence is eliminated 
only after we rescale all variables by A 2 . The implication 
is that the relevant f , r, and 4> are all of the order of A 2 if 
the dimensionless g and u are of order one. Thus, unlike 
the classical case, where the dependence on A was fully 
absorbed into fo, here the low-energy behavior becomes 
non-universal. To proceed - and to make later compar- 
ison with the RG results - here we consider the case of 
fo much smaller than A. This is achieved by taking the 
dimensionless couplings g and it to be small, which also 
implies that r and </> are small compared to A 2 . 

In the absence of the preemptive Ising-nematic in- 
stability, the magnetic instability would again occur at 
fo = 0. Searching for the solution with </> =^ 0, we in- 
troduce as before z — cj)/r, solve the first equation for 
r(z) and substitute the result into the second equation. 
This gives the relation between fo and z. For small g this 
relation takes the form 



r = A 2 e" 
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where 
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Note that fo <C A 2 for g <C 1, as we anticipated. The 
dependence on z comes through /^(z). This function, 
which we plot in Fig. 16, monotonically increases with 
increasing z, implying that the first instability occurs at 
2 = 1. Because f(z) does not depend on a, this result 
obviously holds for any a. 

When g and u are not small, the relation between Tq 
and z is more complex, but the result is the same - the 



first instability upon decreasing fo is into the state with 
the largest <f> = r. As in earlier cases, an instability with 
4> = r implies that the Ising and magnetic transitions 
occur simultaneously and that both are first-order tran- 
sitions. We analyzed the whole phase diagram and again 
found that there is only one transition line at which both 
<j) and M jump to finite values. At smaller fo, <fi an d M 
monotonically increase. 



2. d = 3 

For d = 3, the dependence on the upper cutoff is more 
severe - it is power-law rather than logarithmic. The set 
of two equations on r and cj> becomes, after we integrate 
over frequencies 
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The dependence on A can be eliminated by rescaling 
(r,cp,f Q ) -> A 2 (r,0,f o ) and (g,u) -> (g,u)/(2%A). In- 
troducing, as before z = <f)/r, we rewrite Eqs. (59) in 
rescaled variables as 



q\pr / du^/ut&nh 1 
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(60) 



One can easily make sure that the the first equation in 
(60) has a solution only when the rescaled g is above the 
threshold g cr = 1/2. Once g is above g cri the relevant 
rescaled fo is of order one, i.e. the actual fo is of or- 
der A 2 . While the model remains perfectly well defined, 
universal predictions with respect to the low energy be- 
havior cannot be made. Still, like in d = 2, one can make 
the relevant fo to be much smaller than A 2 . For this, one 
has to place g close to the critical value, g = 1/2 + e, and 
consider e <C 1/A. Expanding in e in (60) and relating r 
to z we obtain, in the original variables, 



r = -A A {g - g cr ) (1 + a) f 3 (z) 



(61) 



where 
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This function, which we plot in Fig. 17, increases with 
z such that the first instability occurs into the state with 
z = 1, i.e. cf> = r. This implies that the Ising-nematic and 
magnetic orders appear simultaneously, via a first-order 
transition. We see therefore that at T = there is no 
difference between d — 3 and d = 2 - a first-order simul- 
taneous Ising/magnetic transition occurs in both cases. 
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Figure 17: (Color online) The function f$(z) from Eq.(62). 

3. 2 < d < 3 

We also analyzed the case of arbitrary d between 2 
and 3 and found the same result as for d = 2 and d = 
3: there is a first-order Ising-nematic transition which 
triggers a simultaneous first-order magnetic transition. 
One can easily show that this result holds for all d > 1 
in the quantum limit, and for all d > 3 in the classical 
limit. The analogy between quantum systems in d > 1 
and classical systems in d > 3 is not surprising because 
the dynamical exponent is Zd yn — 2, meaning that the 
behavior of a quantum system in d dimensions is the 
same as that of a classical system with d + 2 dimensions. 



4- Anisotropic d = 3 case 

The same result - a simultaneous first-order nematic 
and magnetic transition, holds also for anisotropic d = 3 
systems with the susceptibility \q given by Eq. (47). In 
this respect, at T = there is no difference between the 
system behavior in 2 < d < 3 and for anisotropic 3D 
dispersion, no matter what the degree of anisotropy is. 
We stress again that this conclusion is not universally 
true and that it could be possible to construct models 
with same order parameter symmetry that display second 
order transitions (see, for instance, Ref. [57]). 



C. Phase diagram at arbitrary T 

We now combine the quantum and classical analysis 
and consider the phase diagram at a finite T7 ~ fo when 
both quantum and classical fluctuations are equally rel- 
evant. 



1. d = 2 

For d — 2, we have at high temperatures a second 
order Ising transition at a > 2 and a first-order Ising 
transition into 4> < r at 1 < a < 2. At T = 0, we have 



instead a first order transition into <j> = r and rh 7^ for 
all a. A simple analysis shows that the behavior at any 
finite T remains the same as at high temperatures simply 
because at any T > classical fluctuations do not allow 
a non-zero magnetic order. The value of a at which the 
first-order Ising transition becomes second order changes 
with T, but the phase diagram at any T still consists 
of a single line along which the system undergoes either 
first-order or second-order Ising transition. 



2. d= 3 

At d = 3, the behavior at large T and at T — are 
identical - in both cases the first instability is into a state 
with (j) — r (a simultaneous first-order Ising/magnetic 
transition). A simple analysis shows that this behavior 
holds for any T, no matter how small or large. 



3. 2<d<3 

This case is the most interesting one. At high temper- 
atures, all three types of transitions are realized, depend- 
ing on a, while at T = the system only undergoes a 
phase transition into a state with cj> = r. As a result, at a 
given d between 2 and 3, the character of the transition 
changes as a function of a at a fixed T, and as a function 
of T at a fixed a. 

We verified that the phase diagrams do not change if we 
impose upper cutoff on the frequency summation rather 
than on the integration over momentum. The former is 
more convenient for numerical calculations, and below 
we use frequency rather than momentum cutoff. Phys- 
ically, the frequency cutoff A„ becomes more important 
than the momentum cutoff A = A q if the frequency de- 
pendence of the bosonic x(q, v n ) becomes stronger than 
7|f„| at energies smaller than A q . 

For a generic 2 < d < 3 and an arbitrary T, the equa- 
tions for r and <j) become, after integrating over momen- 
tum 
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(63) 



where u and g are defined in (37). The frequency summa- 
tion extends up to n = n max = A u /(2nT). Introducing 
as before z = (f>/r and rescaling in addition the temper- 
ature T — "fT/f , we numerically extract r as a function 
of z from the second equation, substitute into the first 
equation, and obtain fo as a function of z. In Fig. 18 we 
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Figure 18: (Color online) The phase diagram as a function of 
the ratio of the couplings a at fixed reduced temperature T 
and dimensionality d (panels (a)-(b)), as a function of d at 
fixed a and T (panels (c)-(d)), and as a function of T at fixed 
a and d (panels (e)-(f)). The parameters a and T are defined 
in the text. We set n maa; = 100. 



plot fo(z) and the phase diagrams upon varying a, d and 
T. We see the same trend as in the classical phase dia- 
gram, namely, as a gets larger at some fixed T and d, the 
Ising and magnetic transitions split, and at even larger a 
the Ising-nematic transition eventually becomes second 
order (panels (a)-(b)). The same trend holds upon the 
reduction of dimensionality (i.e., make the system more 
two-dimensional) at a fixed a and T (panels (c)-(d)) and 
also upon decreasing T at a fixed a and d (panels (e)-(f)). 

Overall, we see that the phase diagram as a function 
of a does not change qualitatively between high temper- 
atures, when the classical approximation is valid, and 
T = 0(1), when both classical and quantum fluctuations 
are relevant. The only real difference is the change in the 
critical values a c \ at which the Ising-nematic and mag- 
netic transitions split, and a c i at which the Ising-nematic 
transition becomes second-order. As T gets smaller, both 
critical a become larger. Note that the rescaling of T to 



T = "fT/fo is important for obtaining the correct tem- 
perature dependence of the transitions due to hidden T 
dependence in z = <f>/r via r = r(T). 

4- Anisotropic d = 3 case 

We analyzed the two cases of anisotropic magnetic dis- 
persion, Eqs. (47) and (48), at a finite T. In each case 
we found the same behavior as at large T, with different 
phase diagrams depending on the degree of anisotropy, as 
discussed in Fig. 12. For strongly anisotropic, quasi-2D 
spin susceptibility Xqi the results are similar to the case 
2 < d < 3, i.e the magnetic tri-critical point is located 
on the left of the nematic tri-critical point, and in the 
intermediate region the nematic transition is first-order, 
while the magnetic transition is second-order (recall Fig. 
18). For weaker anisotropy, however, the two tri-critical 
points interchange, as we observed in the purely classi- 
cal analysis, see Fig. 12. Once this happens, in the in- 
termediate region between the two tri-critical points the 
nematic transition is second-order, while the magnetic 
transition is first-order. The only difference with respect 
to the purely classical case is that both a c \ and a C 2 shift 
to larger values at a smaller T. 



IV. RG ANALYSIS AT A FINITE N 

The mean-field analysis is quite straightforward, but it 
is rigorously justified only in the artificial limit of large 
JV, where, we remind, N is the number of components of 
the A fields. The actual number of spin components is 
N = 3, and it is by no means guaranteed that the behav- 
ior at N = 3 is the same as at large N. To verify this, we 
need to return back to the effective action S e g [Ax, Ay], 
Eq. (7), and use a complementary approach which is not 
restricted to large N. One such approach, commonly 
used to study phase transitions, is the renormalization 
group (RG) technique. In RG, one progressively inte- 
grates out contributions from high energies down to E 
and analyzes how the parameters of the effective model 
vary with L = log W/E, where W is the bandwidth (the 
highest energy scale in the problem). In our case, the pa- 
rameters are u and g, and in the RG approach one studies 
the flow of the running couplings mj, and gi, and of any 
other coupling generated by the RG flow. Alternatively, 
one can vary the distance to a transition, i.e. vary fo, in 
which case L = log mW/fo ■ 

Quite generally, the RG flow may lead to three types of 
behavior depending on the structure of the RG equations 
and on the bare values u and g. One possibility is that 
the couplings Ul and gh flow to zero, which implies that 
there is no preemptive Ising-nematic transition. Another 
possibility is that g^ and ul flow to infinity (more accu- 
rately, to strong coupling) in such a way that the stability 
condition for the effective model «/, > g^ is not broken. 
In this situation, the system undergoes a second-order 
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Figure 19: (Color online) The one-loop diagrams responsi- 
ble for the renormalization of the two bosonic vertices un = 
(u — g) /2 and U12 = (u + g) /2, with their respective combi- 
natorial factors. The wavy lines correspond to the interac- 
tions, whereas the continuous lines correspond to Ai and A2 
(red/light gray and blue/dark gray lines, respectively). No- 
tice that the two first diagrams in each line contain a closed 
bosonic loop and, therefore, an overall pre-factor of N. 



Ising-nematic transition at the scale L cr = logmW/fQ r 
at which g^ diverges. The third possibility is that the sta- 
bility condition ul > gL gets broken at some L* < L cr . 
In this case, the effective action becomes unstable with 
respect to a discontinuous variation of A^-y, and the 
system undergoes a first-order transition. 

The RG approach is still a weak-coupling approach in 
the sense that the bare couplings u and g are assumed 
to be much smaller than the bandwidth. The advan- 
tage of the RG technique is that it can be applied to 
any TV, and, from this perspective, it goes beyond mean- 
field approximation. However, the RG approach has its 
own limitations - it can be rigorously justified only in 
the marginal dimension d e ff = d + Zd yn = 4, when the 
renormalizations are logarithmic. In our case Zdyn = 2, 
meaning that the marginal behavior takes place at T = 
and d — 2. Still, we obtained the N = 00 phase diagram 
at T = and d = 2 in the previous section, and it is 
instructive to compare that phase diagram with the RG 
phase diagram at arbitrary N to verify whether or not 
the behavior at ./V = 3 is the same as at N = 00, at 
least in this particular case. We remind that our N = 00 
quantum analysis in d = 2 shows that the system under- 
goes a first-order transition into the magnetic state with 
<j> = r for any a > 1. We now analyze how the phase 
diagram looks like for arbitrary N. 

We derived the one-loop RG equations for the flow of 
the running couplings ul ad g^ by the momentum-shell 
method (see, for instance, Ref. [58]) and also derived 
them by analyzing the parquet diagrams for the renor- 
malization of the four-boson vertices, which are presented 
in Fig. 19 together with their respective combinatorial 
pre-factors. Notice the special role of the diagrams which 
contain a closed bosonic loop. Summation over the in- 
ternal bosonic indices yields a factor N, which does not 
appear in the diagrams without closed loops. The mean- 
field results are reproduced if we consider only these di- 
agrams. 



For further convenience, we rescale the coupling con- 
stants to u, g — > 2N (u, g). Combining logarithmic contri- 
butions from all parquet diagrams we can cast the renor- 
malizations of u and g into the form of differential RG 
equations 
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(64) 



where Xl = dX^/dL. 

For the particular cases N — 3 and N — 6, these equa- 
tions reduce to those obtained in Refs.[34] and [57] for 
the similar model. We checked that no other couplings 
allowed by symmetry are generated by the RG flow if 
only u and g have non-zero initial values. In particular, 
if the third coupling v in Eq. (7) is initially zero, it is not 
generated by RG. 

We recall that both ul and gL are initially positive 
(bare values u and g are given by Eqs. 8), and that the 
bosonic action (7) is stable as long as Ul/ql > 1- If this 
condition breaks down at some L, one of the coefficients 
of the quartic terms becomes negative and the system 
undergoes a first-order transition into a state with a non- 
zero <f>. 

For N — > 00 the equations for ul and gL decouple and 
we can easily solve them and obtain 
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We see that it£ flows to zero, while g^ increases un- 
der the RG flow, diverging at L cr — 7^. If g^ would be 
the only parameter in the problem, this divergence would 
indicate a preemptive second-order Ising-nematic insta- 
bility, since the susceptibility of the Ising-nematic order 
parameter diverges at L cr . However, in our case there 
are two couplings, and the action is stable only as long 
as ul > gL- This condition breaks down at a smaller 
^* = 4so — 4«o ' ^ e ^ ore 9 L diverges. The outcome is that 
for N = 00 and d = 2 the system undergoes a first-order 
Ising-nematic transition at T — 0. This is in agreement 
with the mean-field analysis. 

At a finite N, the two equations are coupled and both 
ul and gL can diverge. To understand what happens in 
this case, it is convenient to define the ratio = ul/ 9l 
and re-express the flow equations in terms of gL and pl'. 
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It is straightforward to verify that this set of RG equa- 
tions has several fixed trajectories along which pl is a 
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constant and gL evolves. The fixed trajectories are ob- 
tained by setting pl — in Eq. (66). Solving the 
quadratic equation (— 1 + 2/N)p 2 L — pl — 2/N = we 
find two fixed values p\ = —1 and p^ = —2/ (N — 2). 
The coupling gi diverges along the fixed trajectory with 
PL = pi as <?l = 2g^ (1 + 8 / N) . Along the sec- 
ond fixed trajectory p2, gh evolves according to g^ = 
2g\ (N 2 + 8) / (N 2 - 2N) , i.e. it diverges for TV > 2 and 
tends to zero for N < 2. 

To understand which trajectory is stable and which is 
not, we consider small deviations from a fixed trajectory, 
Pl = Pi + SpiL, and expand the flow equations to lowest 
order in SpiL- We obtain 



SpiL = 2g L SpiL 

Sp2L = ~2gLSp2L 
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(68) 



We see that p\ is a stable fixed trajectory for N < 
4 and unstable for N > 4, while p2 is a stable fixed 
trajectory for N > 4 and unstable for N < 4. 

There is also the third fixed trajectory gL = ( p% — 
oo). Expanding near gL — we obtain from Eq. (64) 
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whose solution is 
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The fixed trajectory gL = is stable as long as gL 
remains small compared to ul- Evaluating pl — UL/gL 
from the solutions of (70) we find 
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We see that, for N > 2, pl decreases under the RG 
flow such that eventually gL exceeds ul- This implies 
that the fixed trajectory gL — is unstable. For N < 2, 
however, pl increases under the RG, and the trajectory 
gL = is stable. 

Combining this analysis with the numerical solution of 
Eqs. (66) and (67) at intermediate energies, we obtain 
three different regimes of system behavior depending on 
the value of N. 

For N > 4, the RG trajectory is as shown in Fig. 20 
(a)-(b). For arbitrary a — u/g, the system approaches 
the stable fixed trajectory Pl = P2 = —2/ (N — 2). Since 
P2 < 0, and the bare value of p is positive and larger 
than one, the running coupling pl necessarily becomes 




(a) 



(b) 



0.2 






N = 3 








N L 




820 840 






L 

1 1 ^ 

5000 10000 



Figure 20: (Color online) RG flow of the running coupling 
constants u (black line) and g (red/light gray line), as well 
as their ratio u/g (blue/dark gray line), as function of L. 
In (a)-(b), we present the results for N — 6; in (c)-(d), we 
display the results for iV = 3; in (e)-(f) the resufts shown 
are for N = 1. In all cases, the initial values were u = 1, 
g = 0.1. The dashed lines in (b) and (d) refer to the stabfe 
fixed trajectories p2 = —2/ (N — 2) andpi = — 1, respectiveiy. 



one at some L = L* along the RG flow. At this L* 
the action becomes unstable and the system undergoes a 
first-order transition into a state with a non-zero Ising- 
nematic order parameter. The only difference between 
finite N > 4 and N — 00 is that, for a finite N, ul 
eventually flows to infinity while for N — > 00 it flows to 
zero. This difference, however, does not play any role 
in our consideration as the RG flow makes sense only as 
long as UL/gL remains larger than one. 

For 2 < N < 4, the fixed trajectory P2 becomes un- 
stable and cannot be reached if the RG flow starts with 
u > g > 0. The stable trajectory to which the sys- 
tem flows is now pi = —1, as shown in Fig. 20 (c)-(d). 
Near this fixed trajectory we find from Eqs. (68) and 
(67) that gL increases and diverges at some L cr , while 
PL + 1 = 5piL oc {l/g L ) (i - N) ' {N+8) -► 0. 

While the fixed trajectory is different for 2 < N < 4 
and for N > 4, the behavior relevant to our purposes 
remains the same for all N > 2, namely in the process 
of RG flow towards a fixed trajectory, the ratio UL/gL 
reaches one at some L — L* . At this point, the action 
becomes unbounded and the system undergoes a first- 
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Figure 21: (Color online) Structure of the RG flow in the (g, u) 
plane for the representative values (a) N = 6, (b) TV = 3, and 
(c) N = 1. The fixed trajectories u — ~Pig are shown as 
dashed lines in the lower-half planes (pi = — 1 for (b) and (c) 
and p2 — —2/ (TV — 2) in (a)), and are approached only at 
very large L. In the upper-half planes, the dotted line u — g 
separates the second-order and first-order regimes. In (c), 
the dotted line u = 2g/ (2 — TV) separates initial conditions 
that fall in the basis of attraction of the g — fixed point 
(p3 = oo, green/light solid line) from those that flow to the 
pi = — 1 fixed trajectory (red/dark solid line). 



order phase transition into a state with a finite Ising- 
nematic order parameter. The dependence on TV is then 
only quantitative - as TV decreases, the scale L* gets pro- 
gressively larger. To see this, we note that for any finite 
TV, the couplings ul and first decrease under the RG 
flow and only at larger L reverse the trend and approach 
the fixed trajectory at —¥ oo and ul — > — oo. To get 
an estimate of how L* evolves with TV, we use the ap- 
proximate Eq.(71) and identify L* with the RG scale at 
which pl becomes one. We obtain 



L* 



■x 



(72) 



For a given u/g 3> 1, the value L* rapidly increases 
with decreasing TV; for instance, for u/g = 10, the ratio 
of L* for TV = 3 and for TV = 5 is 10 4 . This means that 
at smaller TV the system behaves over a wide range of 
energies as if magnetic fluctuations were absent, and only 
very near T^,o it recognizes that it actually undergoes a 
first-order nematic transition. 

For TV < 2, a new behavior becomes possible, as 
shown in Fig. 20 (e)-(f). Now the fixed trajectory 
P2 = 2/ (2 — TV) crosses the region u > g > from 
where the RG flow begins. Once the bare u and g are 
such that u/g > 2/ (2 — TV), the RG flow is sandwiched 
between the fixed trajectories p2 and P3 which are re- 
spectively unstable and stable for TV < 2. The RG flow 
then moves both g^ and ul towards gL,UL = 0, keeping 
gi < ul, i.e. without crossing the first-order instabil- 
ity line. In this situation, no preemptive nematic insta- 
bility develops, and the system only undergoes a mean- 
field magnetic transition at T)v,o- If, however, the ini- 
tial u/g < 2/ (2 — TV), the system behavior is the same 
as before, with the couplings evolving towards the fixed 
trajectory p\ = —1, and pl becoming equal to one at 
some scale L*, at which the system undergoes a first- 
order nematic transition. In the formal limit TV = 0, the 
whole region u > g falls into the basis of attraction of the 
9l = ul = fixed point, i.e. there is no preemptive ne- 
matic instability for any u/g > 1. The general structure 
of the RG flow equations in the (g, u) plane is shown in 
Fig. 21. 

For d = 2 and T = this scenario of no preemptive 
nematic instability holds only for TV < 2 and does not 
affect our actual case of TV = 3, for which the system 
behavior under the RG flow is qualitatively the same as in 
the mean- field, TV = oo analysis. What happens for d > 2 
and/or a finite T is unclear because, in the absence of the 
logarithmic terms, the approximations leading to the RG 
equations are not justified. It is possible in principle that 
a preemptive nematic transition does not occur for some 
large enough a — u/g. If this is the case, then there must 
be a reentrant behavior for large a in Fig. 8, i.e. the 
magnetic transition temperature Tjv must reverse trend 
and come closer to the nematic instability. However, the 
more likely scenario is that the phase diagram which we 
obtained in the mean-field approximation survives for the 
actual TV — 3 component bosonic field for all a. 
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V. CONSEQUENCES OF THE ISING-NEMATIC 
ORDER 

A. Orbital Order 

ARPES measurements on detwinned samples have 
found that the onset of resistivity anisotropy is accom- 
panied by the onset of orbital order in the paramagnetic 
phase, with different occupations for the d xz and d yz or- 
bitals [10]. One possibility, explored by several authors 
in different contexts, is that this orbital ordering is an 
intrinsic instability of the system [18-24, 26, 27]. In line 
with the theme of this work, we explore another pos- 
sibility, namely that the orbital order is induced by the 
Ising-nematic order. This scenario is generally consistent 
with the small value of the measured orbital polarization. 
To investigate this scenario quantitatively, we consider a 
simplified two-orbital model in which the entire X pocket 
has dyz character, while the entire Y pocket has d xz char- 
acter [30], and assume that there is a splitting A or b be- 
tween the on-site energies of the d xz and d yz orbitals. In 
the presence of such splitting, the fermionic dispersion 
becomes anisotropic, see Fig. 22, and the Hamiltonian 
acquires the additional terms 

"Horb = ~ Aorb ( C kk Q C X,ka - 4- ikQ Cy :kQ ) + y A^ rb 

(73) 

Including these two terms into the Hubbard- 
Stratonovich procedure and expanding the effective ac- 
tion in powers of A or b we obtain 

S[A X , Ay,A orb ] = S eS [A x , A Y ] + ^A 2 mh 

-w (A 2 X - A 2 Y ) A orb (74) 
with S Q g [Ax, Ay] given by Eq. (7) and: 

a = a a + 4 / & Xk 

Jk 

w = 2 f G T , k G\ >k (75) 

Jk 

Evaluating the integrals and expanding around perfect 
nesting we find a > and w = —cfi with c > 0. The w 
term describes the linear coupling between the orbital 
and Ising-nematic order parameters, i.e. the develop- 
ment of one order triggers the development of the other. 
Differentiating Eq. (74) with respect to A or b we obtain 

(A olh ) = ^(A%-A^) (76) 

Since w cx —fi and /x scales with doping, the in- 
duced orbital order is expected to be small at small dop- 
ing. This is in accordance with the experimental data. 
The sign of the orbital splitting also agrees with the 




Figure 22: (Color online) Schematic representation of the ef- 
fect of orbital ordering on the Fermi surface (left panel) and on 
the band dispersions (right panel). The light /orange curves 
have d xz orbital character, while the dark/blue curves have 
dy Z orbital character. The solid lines correspond to the Fermi 
surface and band dispersions without orbital order (T > T s ), 
whereas the dashed lines refer to temperatures below the on- 
set of orbital order (T < T s ). In this figure we considered 
A orb < 0. 

data[10]: polarized ARPES measurements on detwinned 
electron-doped Ba (Fei_ 2; Co 2: ) 2 As2 samples have shown 
that when (A^-) > (Ay), (A or b) is negative. This is 
consistent with (76) since for electron-doped materials 
ji > and w < 0. To the best of our knowledge, similar 
measurements have not been carried out for hole-doped 
samples. Interestingly, Eq. (76) shows that the sign of 
(Aorb) should change for hole-doped materials, which are 
described in our model by ji < 0. Note that long-range 
magnetic order can also induce orbital polarization [22]. 

B. Structural order 

The same reasoning also applies to the interplay be- 
tween the Ising-nematic and the orthorhombic orders. 
The structural order is detected experimentally as the 
difference between the lattice constants a and b along the 
x and y directions of the Fe-plane, respectively. In the 
ideal situation, structural order appears only below a par- 
ticular structural transition temperature T s . In reality, 
however, some orthorhombic distortion can be present at 
any T due to internal strains. In the case of detwinned 
samples, a small in-plane strain is applied along one of 
the orientations [4, 6, 59, 60]. Then, strictly speaking, 
e s = a — b is never zero, i.e. there is no well-defined T s 
for finite strain. Still, experimentally one can identify the 
crossover temperature below which e s sharply increases. 

The relationship between Ising-nematic and structural 
orders can be obtained in the same way as in the previous 
subsection. Introducing the orthorhombic order param- 
eter in a way similar to Eq. (73) yields 

k 

(77) 

where A is a coupling constant and C s is the shear 
modulus. Including these two terms into the Hubbard- 
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Stratonovich procedure we obtain, in the mean-field ap- 
proximation, the same effective action as Eq. (74), but 
with e s instead of A or b ail d with the renormalized shear 
modulus C s instead of a. Accordingly, the orthorhombic 
and Ising-nematic order parameters are linearly related: 

(e.) = ^(A* x -A$r) (78) 

i.e. one order immediately triggers the other. Notice, 
however, that in distinction to the case of orbital or- 
der (Eq. 76), the proportionality constant between (e s ) 
and the nematic order parameter contains not only the 
parameter w, but also the magneto-elastic coupling A, 
which can itself depend on the chemical potential and on 
additional details of the band structure. 

The linear relation between the two orders has been 
discussed in a number of earlier papers [32, 35, 61], and 
is not surprising because both orders break the same Zi 
symmetry between the x and y directions (the orbital or- 
der does the same). It also implies that, in detwinned 
samples, the applied strain plays the role of a conjugate 
field to the Ising-nematic order parameter [40]. We recall 
that the proportionality coefficient w scales with /i and 
is generally small. As a result, if the orthorhombic or- 
der is induced by the Ising-nematic order at T = T s , the 
orthorhombic order parameter is initially small and may 
become visible only at some distance below the Ising- 
nematic transition. This may explain why recent mag- 
netic torque experiments observed nematic order up to 
higher temperatures than the structural order [11]. We 
emphasize, however, that there is only one well-defined 
transition temperature below which the tetragonal sym- 
metry is broken. 

Note that the coupling to structural degrees of freedom 
also renormalizes the nematic coupling constant as [35] 

9 = 9 + — (79) 

Therefore, even if initially g = 0, the coupling to the 
lattice generates a non-zero nematic coupling and, con- 
sequently, nematic order [61]. 

C. Pseudogap behavior 

Recent ARPES data on NaFeAs, whose strongly 
anisotropic Fermi surface is very similar to that of the 
1111 compounds [62], show that the reconstruction of 
the fermionic spectrum begins at temperatures around 
T s rs 54 K, where the nematic order sets in, rather than 
at T/v ~ 39 K, where the stripe magnetic order develops 
[36]. Such a reconstruction increases below Tn, which is 
a good indication that the effect likely has a magnetic 
origin. 

Thermal magnetic fluctuations do give rise to a pseu- 
dogap behavior (often termed as magnetic precursors) by 
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Figure 23: (Color online) Inverse squared magnetic correla- 
tion length £ -2 oc r (in arbitrary units), as function of the 
reduced temperature Afo oc T — T s , for the case a = 4 in the 
phase diagram of Fig. 8. The dashed lines denote the split 
second-order nematic (T s ) and magnetic (Tn) transitions. 
Note the sharp increase and the discontinuity in the deriva- 
tive of £ at the nematic transition. The insets show schemat- 
ically, at different temperatures, the frequency-dependence of 
the fermionic spectral function Ahs (w) at the hot spots, i.e. 
the Bragg points of the band dispersion er,k = £(x,Y),k+Qi ■ 
Notice the development of peaks at cj — ±Asdw below the 
nematic transition, where Asdw denotes the value of gap that 
is opened due to SDW order at zero temperature. 



transferring spectral weight from small frequencies to a 
frequency comparable to the spectral gap developed in 
the magnetically ordered state [63]. The only precondi- 
tion for magnetic precursors is that the magnetic corre- 
lation length £ must be large enough. 

How can the pseudogap develop at the nematic tran- 
sition? At first glance, this seems a mere coincidence 
because the nematic order is a collective instability in 
the particle-hole channel with momentum q = 0, and a 
q = boson cannot reconstruct the fermionic spectrum. 
However, one of our key results is that at the nematic 
transition the magnetic correlation length increases ei- 
ther discontinuously or very sharply. Even away from the 
nematic tri-critical point, where the Ising-nematic transi- 
tion is second-order, the magnetic correlation length has 
a discontinuous derivative at T s , since £ -2 — > — 4>, 
making it increase faster. This behavior is shown in Fig. 
23, where we plot the inverse magnetic correlation length 
r oc as a function of temperature for a — 4 for the 
phase diagram of Fig. 8. Interestingly, in NaFeAs, as 
well as in other iron pnictides [64] , a significant enhance- 
ment of magnetic fluctuations was observed [65, 66] just 
below T s by nuclear magnetic resonance (NMR). Once 
the correlation length jumps (or sharply increases) to a 
larger value, the strength of thermal magnetic fluctua- 
tions rapidly increases. As a result, the fermionic spectral 
function A\ (to) develops a magnetic pseudogap via the 
transfer of spectral weight from zero to finite frequencies 
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at the hot spots er,k = £(x,Y),k+Qi ( see Fig. 23), what 
leads to the reconstruction of the fermionic spectrum, al- 
though zero- frequency states appear only below Tjv . We 
argue therefore that the jump (or sharp increase) in £ at 
the nematic transition is the "glue" that links the nematic 
instability and the development of the pseudogap in the 
fermionic spectral function. 



VI. COMPARISON WITH EXPERIMENTS AND 
EARLIER THEORIES 

A. Experimental phase diagrams 

We now compare our results of Section III for 
anisotropic systems to the experimental phase diagrams 
of the 122 and 1111 compounds. Our theoretical phase 
diagram is plotted in Fig. 1 (a) and (b). Our goal is to 
relate the changes of the system behavior as function of 
the parameter a = u/g to the measured changes imposed 
by doping, pressure, and alkaline-earth substitution. For 
this we need to (i) properly place the parent compounds 
onto our phase diagram, (ii) decide which of the two sim- 
ilar but not identical phase diagrams in Fig. 1 is more 
appropriate, and (iii) understand how a varies with the 
experimental parameters. We first consider 122 materials 
and then briefly discuss 1111 materials. 

We begin with the parent 122 materials. There 
are three well studied types of 122 systems, namely, 
CaFe 2 As 2 , SrFe 2 As2, and BaFe 2 As2- The first two mate- 
rials undergo a strong simultaneous magneto-structural 
first-order transition, as evidenced from several thermo- 
dynamic measurements [67-69]. This clearly places both 
materials well into the region of small a, to the left of 
the first tri-critical point - either a c \ in Fig. 1(a) or a C 2 
in Fig. 1(b). For BaFe2As2, X-ray diffraction [15, 70] as 
well as high- accuracy magnetization measurements [16] 
find a very small splitting between the structural and 
magnetic transitions (T s ~ 141 K, T N ~ 140 K). This 
places BaFe2As2 in the region of larger a (i.e. smaller 
g) , which can be either slightly to the right of the mag- 
netic tri-critical point a cl in Fig. 1(a), or slightly to 
the right of the nematic tri-critical point a C 2 in Fig. 
1(b). A smaller a for the parent compounds CaFe2As2, 
SrFe 2 As 2 compared to the parent compound BaFe 2 As 2 
can be explained by the difference in the values of the 
shear modulus C s . The analysis of their phase diagrams 
under pressure [71-73], as well as of their mechanical 
properties [74], shows that CaFe2As2 and SrFe2As2 are 
softer than BaFe 2 As 2 . According to our Eq. (79), softer 
systems with smaller shear modulus C s (and possibly 
larger magneto-elastic coupling A) have larger g, and 
hence smaller a. 

Which of the two phase diagrams in Fig. 1 is more 
appropriate for BaFe 2 As 2 is a more subtle issue. The 
x-ray data [15, 70] show that the orthorhombic order pa- 
rameter evolves continuously immediately below T s and 
then jumps at the same temperature where the magnetic 



transition takes place, as evidenced by the magnetization 
data [16]. Neutron diffraction data [75] do not detect crit- 
ical magnetic fluctuations above Tat, consistent with the 
idea that the magnetic transition is first-order and simul- 
taneous to the meta-nematic transition. This favors the 
phase diagram of Fig. 1(b). However, the same neutron 
data do not detect a clear jump of the magnetic order pa- 
rameter at the magnetic transition, as expected at a first- 
order transition. So, it is possible in principle that the 
continuous evolution of the nematic/orthorhombic order 
parameter immediately below T s is a secondary effect - 
due to internal strain, for instance. Then, the true ne- 
matic transition would be first-order, while the magnetic 
transition occurring at 7jv < T s would be second order. 
This would be consistent with the phase diagram of Fig. 
1(a). In any case, however, the key observation is that 
the parent compound BaFe2As2 is in the region where 
the structural and magnetic transitions are quite close to 
each other and at least one of them is first order. 

To investigate the effect of doping and pressure, we 
need to understand how a — u/g changes with the chem- 
ical potential \i and the mass anisotropy of the electron 
pockets 5m — m (m x — m y ) /2m x m y . We use Eqs. (8) 
and evaluate u and g for small chemical potential [i/eo 
and small 5m. For the ratio u/g we obtain 




Electron doping adds carriers to the electron pockets 
and increases the magnitude of the chemical potential \x. 
This, according to (80), increases a, i.e., under electron 
doping the system should move to the regime of split 
second-order transitions (see Fig. 1). This is the main 
prediction of our theory. 

This prediction generally agrees with the experimental 
phase diagrams of electron-doped 122 materials, where 
Fe is substituted by Co, Ni, Cu, Pd, or Rh [76]. In 
particular, for Ba (Fei_ 2; Co 2 ;) 2 As2, magnetic and x-ray 
measurements demonstrated [15, 16] that the magnetic 
and structural transitions rapidly split and both become 
second-order above x w 0.022. This is particularly ev- 
idenced by the fact that the peak in the derivative of 
the magnetic susceptibility \ (indicative of a first order 
magnetic transition) is strongly suppressed beyond this 
doping value [16]. The splitting of both transitions upon 
doping has also been observed in Ca (Fei_ x Co x ) 2 As 2 at 
x w 0.039 by neutron diffraction data [17]. In these sys- 
tems, the structural transition remains first-order at least 
for some doping range after the splitting, what is consis- 
tent with the phase diagram in Fig. 1(a) if we iden- 
tify x w 0.039 with a c \. The behavior of the isovalently 
doped compound Ba (Fei_ x Ru x ) 2 As 2 is also consistent 
with our theory. In this material, the chemical potential 
does not change with doping [77], hence the structural 
and magnetic transitions should remain very close for all 
x. Thermodynamic measurements did indeed find that 
the two transitions do not split upon increasing x [78]. 
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The interplay between structural and magnetic tran- 
sitions in the hole-doped (Bai-^K^) Fe2As2 compounds 
is not so well established, with conflicting reports of 
either simultaneous first-order transitions [79], or split 
first-order structural transition and second-order mag- 
netic transition [80]. In the context of our model, hole 
doping adds an extra complication since at least some 
effects of doping are absorbed into the changes of an ad- 
ditional hole pocket at (tt,tt) [81]. This suggests that the 
chemical potential in our effective 4-band model changes 
at a slower rate than in the case of electron doping, since 
extra holes do not necessarily go to the central pockets. 
If this is the case, then a first-order structural transition 
extends to larger dopings, in agreement with the data. 



Consider now the evolution of a with pressure. Band 
structure calculations show that pressure reduces the 
nesting features of the Fermi surface [82, 83], what in 
our model implies that 8m increases. According to Eq. 
(80), a oc l/(Sm) 2 then decreases, bringing the system 
deeper into the regime of simultaneous first-order struc- 
tural and magnetic transitions (see Fig. 1). This agrees 
with the experimental results that structural and mag- 
netic transitions do not split under pressure in ^4eFe2As2 
(Ae = Ba, Co, Sr) [67, 73]. 



We now briefly consider the 1111 materials, whose 
Fermi surfaces are also similar to the previously discussed 
NaFeAs compound. In these systems, the nematic and 
magnetic transitions are split and second-order already in 
the parent compounds. In our theory, these compounds 
should then be placed to the right of a C 2 in Fig. 1(a) 
or to the right of a c \ in Fig. 1(b). The reason why a 
is larger in the 1111 materials is the significantly larger 
degree of out-of-plane anisotropy in the 1111 materials 
compared to the 122 materials, what brings the former 
closer to the d = 2 limit. In our modeling, these systems 
are then described by a smaller effective d. According to 
our theory, as d gets smaller, both tri-critical points ot c \,2 
shift towards smaller values (see Fig. 18 and Eq. (40)), 
extending the regime where the magnetic and structural 
transitions are split and second-order. A similar argu- 
ment was given in Ref. [34]. 



Finally, for the iron chalcogenides FeTei_ x Se x , our 
model is applicable in the regime of intermediary Se dop- 
ing, near the superconducting dome of the (%, T) phase 
diagram. In this region, ARPES measurements [84] re- 
veal that the electronic structure is similar to the one 
considered in our model (see Fig. 2), and neutron scat- 
tering shows [85] that magnetic fluctuations are peaked 
at Qi = (-7T, 0) and Q2 = (0,7r). On the other hand, 
our model is not suitable for the undoped FeTe sample, 
where the same ARPES data show the absence of elec- 
tron pockets centered at Qi and Q2 [84]. 
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Figure 24: (Color online) Nematic order parameter cf> 
(square/blue symbols) and re-scaled magnetic order param- 
eter M (triangles/red symbols) for d = 2.5 and a — 1.8 as 
function of — Afo = ro, cr — ro oc (T s — T) (the same as in Fig. 
9d). Here, M = 0.87M + 0.33. 



B. Linear relation between the magnetic and 
nematic order parameters 

Several experimental groups showed that in some 122 
materials, most notably SrFe2As2, the magnetic and 
nematic order parameters have very similar temper- 
ature dependencies below the simultaneous first-order 
magneto- nematic transition [79, 86, 87]. In our analysis, 
the relationship between M and tfi is given by Eq. (41) 
and can be readily seen in Fig. 9d for a simultaneous 
first-order transition. Re-plotting in Fig. 24 the nematic 
order parameter (f> and the properly re-scaled magnetic 
order parameter M as functions of temperature, we see 
that both follow the same trend, indicating that the re- 
lationship between the two order parameters is nearly 
linear, in agreement with the experimental data. Even 
better linear relation is obtained for a closer to one, as 
we show in Fig. 25, where </>/A</> is plotted explicitly as 
function of M / AM (A0 and AM denote the values of 
the jumps across the first-order magneto-nematic transi- 
tion) . 

At first glance, this near-linear relation seems non- 
trivial, because for split second-order magnetic and ne- 
matic transitions, a straightforward expansion leads to 
4> — cf>(m = 0) oc M 2 . It can, however, be easily un- 
derstood analytically by expanding Eq. (41) around the 
jumps Acj) and AM. We obtain 



A(f> 



A(j) 



(M 



V AM 



AA/V 



(81) 



For small <fi — Ac6 and M — AM, the relationship is 
indeed linear, since a\ is non-zero for all 2 < d < 3 
(see Appendix D). In fact, in SrFe2As2, where the linear 
relation was experimentally observed[87], the measured 
temperature dependent orthorhombic and magnetic or- 
der parameters are rather small compared to the magni- 
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Figure 25: (Color online) Nematic order parameter cj> as func- 
tion of the magnetic order parameter M for d — 2.5 and 
a — 1.05. Note the near linear behavior between them. Here 
Acj> and AM refer to the jumps at each transition. 



tilde of the jumps. Interestingly, when a — > 1 and the 
first-order transition gets stronger, the coefficients of Eq. 
(81) satisfy ai/a 2 — 2 and Oj +2 = for any 2 < d < 3, 
i.e., the dependence has only linear and quadratic terms. 
In the same limit, A(j> and AM become large, since both 
scale as l/(a — 1). Thus, a plot of <f> versus M in absolute 
units would show a strictly linear dependence for a — > 1 . 



C. Ji — J2 and phenomenological models 

The possibility of an Ising-nematic order induced by 
magnetic fluctuations was first proposed for the iron 
pnictides within the localized-moment scenario [31-33], 
built upon earlier results by Chandra, Coleman, and 
Larkin [28]. In this localized-moment approach, one 
considers spins on a square lattice interacting via a 
nearest-neighbor antiferromagnetic exchange J± and a 
next-nearest neighbor antiferromagnetic exchange J 2 . It 
was argued [31] that the hybridization between the Fe 
and As orbitals in the iron pnictides gives rise to a rather 
large J 2 , which can exceed Ji/2. Once this happens, 
the magnetic ground state develops the stripe order with 
0(3) x Z-2 order parameter manifold. 

As shown in Ref. [28], one can interpret the stripe 
order as composed of two inter-penetrating Neel sublat- 
tices, with staggered magnetization Mi and M2. The 
configuration with Mi parallel (anti-parallel) to M2 cor- 
responds to the (ir, 0) ((0,7r)) state. At the mean-field 
level, the two sublattices are uncoupled. Thermal and 
quantum fluctuations, however, induce a coupling be- 
tween the two sublattices that favors the collinear con- 
figurations. Once this coupling is included, the effective 
action takes the form 

Sj x -j* [Mi] = J 2 fq 2 (Mi,, ■ Mi,_ s + M 2 , g • M 2 ,_ g ) 

J q 

+Ji / q x q y (Mi >9 • M 2 ,_ g ) 

J a 



where the magnitudes of Mi and M 2 are fixed, and £ 
is a dimensionless constant which is non-zero because of 
thermal and/or quantum fluctuations. This constant is 
small in a l/S expansion, where S denotes spin, and 
remains small numerically even for S = 1/2 (see, for 
instance [88]). 

It is straightforward to make connection between our 
itinerant model, Eq. (7), and the J\ — J 2 model. First, 
the relationship between the real-space order parameters 
Mi and M 2 and the momentum-space order parameters 
Ax and Ay is A x = Mi + M 2 and Ay = Mi - M 2 . 
The scalar product Mi • M 2 is then the same as A\ 2 X — 
Ay and the nematic coupling is given by g — C,J\j J^- 
Second, in the itinerant approach, the hard constraint 
Mf = M 2 = 1 is replaced by the quartic terms, which 
play the role of soft constraints. We see that the itinerant 
and the J\ — J 2 models are indeed quite similar, and in 
both models the Ising-nematic order results from the Z 2 
degeneracy of the stripe magnetic ground state. 

Whether the two models have identical phase diagrams 
and show the same behavior upon doping, pressure, and 
alkaline-earth substitution is a more subtle issue. In the 
J\ — J 2 model, the analog of a is J|/(C^i)- Taken at 
face value, this quantity is large (because £ is small) and 
depends only weakly on doping and pressure. For in- 
stance, the parameters used in Ref. [32] yield a ~ 100, 
which places both undoped and doped iron pnictides in 
the regime of split second-order transitions, i.e. region III 
of the phase diagrams of Figs. 8 and 14. On the other 
hand, in the itinerant model, a is generally of order one 
and changes with doping and pressure due to the changes 
in the Fermi surface. 

It is also unclear (chiefly due to the lack of results) 
whether the phase diagram of the J\ — J 2 model con- 
tains an intermediate phase in which one transition is 
first-order and the other is second-order, as in Fig. 1. 
The mean-field analysis of Eq. (82) for a quasi two- 
dimensional system did find split second-order nematic 
and magnetic transitions [32] for small £ and a simultane- 
ous first-order transition[57] for large enough £. A similar 
result was obtained by Monte Carlo simulations of Eq. 
(82) in Refs. [89, 90], which treated £ as a phenomeno- 
logical input parameter. It was also shown by mean-field 
[91], RG [34], and Monte Carlo [92] calculations that in 
three dimensional systems with anisotropic magnetic dis- 
persion the degree of anisotropy tunes the system be- 
tween the regimes of split second-order transitions and 
simultaneous first-order transitions. Making the spin in- 
teraction anisotropic in spin space has the same effect, 
i.e. it gives rise to a transformation from split second- 
order transitions to simultaneous first-order transitions 
[93]. However, in all cases, the intermediate phase either 
was not discussed, or was assumed to be absent [34]. If 
only the regimes of split second-order transitions and si- 
multaneous first-order transitions occur in the J\ — J 2 
model, it will be difficult to explain the different char- 
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acters of the transitions observed experimentally upon 
doping, pressure, or alkaline-earth substitution. Also, to 
the best of our knowledge, there have been no studies 
within the J\ — J 2 model of whether the magnetic cor- 
relation length jumps at the nematic transition. As we 
showed in Section V C, such a jump (or a sharp increase), 
obtained with the itinerant approach, is fundamental to 
explain the pseudogap behavior above Tjv. 

Several research groups, using J\ — J2 based models, 
also put forward various arguments to relate the trans- 
formation from simultaneous to split nematic and mag- 
netic transitions to doping and pressure [34, 94, 95]. 
Ref. [34] suggested that the change from simultane- 
ous first-order to split second-order transitions with dop- 
ing is due to doping-induced change of the spin dynam- 
ics from Landau-overdamped (zd yn = 2) to propagat- 
ing (zdyn = !)• However, ARPES measurements later 
showed that hot spots - and therefore Landau damping 
- are present even at optimal doping [96]. The authors 
of Ref. [95] assumed that doping increases the disorder 
concentration and decreases the coupling to the lattice, 
what leads to a decrease in g. However, recent data on 
isovalent doped pnictides [77, 78] show that, when the 
chemical potential remains unchanged, the character of 
the transitions remain the same as in the parent com- 
pound, even for large doping concentrations and large 
in-plane disorder. Ref. [94] suggested that the phase di- 
agram of the iron pnictides is close to a magnetic quan- 
tum tri-critical point, but did not analyze the character 
of the structural transition. 

There have also been studies of nematic and magnetic 
transitions within the itinerant approach to the iron- 
pnictides. Ref. [15] considered a model very similar to 
the one presented here and argued that the phase dia- 
gram of Fig. 1(b) can be obtained even for quasi- two 
dimensional systems if there is a strong enough coupling 
to anharmonic elastic terms. 

A phase diagram similar to our Fig. 1 (b) , with the in- 
termediate regime, was also obtained in Ref. [56]. These 
authors assumed that the structural and magnetic transi- 
tions occur independent of each other, and that the two 
order parameters are linearly coupled, as in Eq. (78). 
The intermediate regime then emerges in some range of 
parameters, primarily due to the presence of a quartic 
term in the elastic free-energy. The authors of [56] , how- 
ever, did not argue why the intrinsic magnetic and struc- 
tural transitions would occur at about the same temper- 
ature. 

One argument of why this may be the case, without 
invoking nematic degrees of freedom, was presented in 
Ref. [97]. The authors of [97] argued that, under spe- 
cial conditions (which may or may not be satisfied in the 
iron pnictides), a SDW instability and a CDW instability 
in the form of orbital currents occur at almost the same 
temperature, with the orbital-current instability occur- 
ring first. They further argued that orbital-current order 
and magnetic order are orthogonal to each other in mo- 
mentum space and coexist. 



As we already mentioned, some elements of the physics 
that leads to our phase diagrams in Fig. 1 are similar to 
the physics of the J\ — J2 and phenomenological mod- 
els. The key elements that distinguish our study from 
previous works are (i) that the whole phase diagram is 
entirely driven by magnetic degrees of freedom, and (ii) 
that the parameters of the effective bosonic model are 
derived from the original itinerant model of interacting 
fermions, such that the evolution of the system behavior 
with doping, pressure, and alkaline-earth substitution is 
fully described within the model itself. 

Finally, we point out that an RG analysis similar to the 
one presented in Section III was performed in Ref. [34] 
for N = 3 components of the magnetic order parameter, 
and in Ref. [57] for N — 6, which refers to the case of in- 
commensurate magnetic order parameter. Interestingly, 
although both works obtained runaway flows, the former 
associated the first-order instability to the divergence of 
the coupling constants, while the latter pointed out that 
it happens much earlier, when the condition ul < 9l is 
first satisfied. Our reasoning is similar to the one in Ref. 
[57]. 



VII. CONCLUDING REMARKS 

We argued in this work that the development of the 
preemptive nematic order and its interplay with the 
stripe magnetic order can be fully understood within 
an itinerant magnetic scenario for the iron pnictides. 
We obtained (rather than assumed) the dependence 
of the nematic coupling upon doping, pressure, and 
alkaline-earth substitution, what enabled us to consis- 
tently explain the character of the magneto-structural 
transitions for a variety of iron pnictides. We also 
demonstrated how the nematic transition induces or- 
bital and structural order, and triggers the pseudogap 
behavior observed in the paramagnetic phase of these 
materials. 
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Appendix A: Microscopic calculation of the 
Ginzburg-Landau coefficients 

In this Appendix, we show explicitly how to calculate 
the Ginzburg-Landau coefficients of Eq. (7) starting with 
the interacting Hamiltonian W = Ho + "Hint of Eqs. (I) 
and (2). First, we introduce the 6-dimensional creation 
Nambu operator: 



J 

-r,kf 



J 



J 



J 



J 



(Al) 



After introducing the bosonic fields Ax and Ay via 
the Hubbard-Stratonovich transformation and evaluating 
the products of Pauli matrices, we can write the partition 
function as: 



Here, Tr (• • •) refers to sum over momentum, frequency 
and Nambu indices. A series expansion in powers of A| 
then gives: 



S cS [A x ,Ay] = ^Tr(g , k Vf + -TY{g . k vf 



— /( 

u spin J x 



A 2 X + Ay) 



Evaluation of the traces yields: 



(A8) 



S cff [A x , A Y ] = J2 n>,A ? + £ u v A 2 A 2 (A9) 

i i,j 

with the coefficients: 



Z = J dA t cMe~ s ^^ 
with the action written in compact form: 



(A2) 



S[*,A i ] = - / *tg-i^ fc + ^_ f (A 2 X + A 2 Y ) (A3) 

Jk "spin Jx 



Here, Aj = | A^ | and the Green's function Q k 1 is given 



by: 



with the bare term: 



Go,k 



Gtm 
Gx,k 
G Y , k 



and the interacting term: 



(A4) 



(A5) 



ro,i = 



+ 2 / Gr,kGi t k 



Uij = / G r ^Gi^Gj^k 

Jk 



(AI0) 



Due to the 7r/2 rotation symmetry relating the X and 
Y bands, it follows that ro,i = ro,2 = ro and uu = U2i- 
After re-arranging the terms, we obtain: 



SW[Ax,A y ] = r (A 2 x + A 2 Y ) 
+ 



Ui2 + Uii \a x + a y ) 2 



u 12 - Uu 



(A x -A Y )lAU) 



with the coefficients u = u\ 2 + «n and g — u\ 2 — Uu 
given in Eq. (8). 



Appendix B: Effects of the anisotropic momentum 
dispersion 



-A x -A y 
V = I -Ax 

-Ay 



(A6) 



Here, we defined the 2x2 matrices G^fe = G^I and 
Aj = Aj • <r, with identity matrix I, Pauli matrices ct 3 , 
and G~l — iuj n — £, ; k the non-interacting single-particle 
Green's functions. 

It is now straightforward to integrate out the fermions, 
since the action is quadratic in them, and obtain the 
effective magnetic action: 



S off [Ax, Ay] =-THn(l-0 o ,fcV) + — f (Ax + Ay 

M spin Jx 

(A7) 



In this Appendix we briefly show that the inclusion 
of an anisotropic momentum dispersion in the bare dy- 
namic susceptibilities Xi.q do not change our main re- 
sults. After denoting the renormalized susceptibilities by 
Xx\ = Xx\ + V> - 4> and Xx\ = X Y \ + 1> + 4>, with 
X^q = ro + fi.q + 7 \v n \ j we can rewrite the mean-field 
equations (17) as: 



= t; f ixx, q + XY, q ) 



lj q (xx, g 



XY,q) 



(Bl) 



In the main text, we considered the case of an isotropic 
momentum dispersion /j q = q 2 . Most generally, /j jq 
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will have an anisotropic form preserving the tetragonal 
symmetry of the system: 



/», q = ql(l± v) 



ll (1 T »?) 



(B2) 



where — 1 < rj < 1 and the upper (lower) sign refers to 
band X (Y). Indeed, inelastic neutron scattering mea- 
surements hnd this form for the dynamic susceptibility 
in several iron pnictide compounds [98, 99]. 

Substituting it in the self-consistent equations (Bl), 
we can re-scale the momentum by q x — q x \/l ± rj and 
% = "F V depending on whether the integral in- 

volves Xx,q (upper signs) or XY.q (lower signs). In ei- 
ther case, the Jacobian of the transformation is the same, 
yielding: 



1> 



rp + tp + q +j\i>m\ 

V 1 - v 2 h ( r o + i> + q 2 + iWm\) 2 - 4> 2 

9 



yl - V 2 h (ro + i> + q 2 + l\v n 



(B3) 



Thus, comparing the previous equations with the orig- 
inal mean- field expressions (17), we conclude that the 
only effect of the anisotropic dispersion is to renormalize 
the coupling constants u and g in the same way, yielding: 



9 = 



9 



(B4) 



This does not change the value of the ratio a — u/g — 
u/g, implying that the phase diagrams discussed in the 
main text remain valid, with the same values for the tri- 
critical points a c i and a c i- The only modification will 
be in the absolute value of the temperature, since Fo is 
proportional to g (see Eq. 38). 



Appendix C: Anisotropic 3D model with quadratic 
dispersion 

In this Appendix we consider the behavior of the mag- 
netic and nematic tri-critical points in the anisotropic 3D 
model with the bosonic susceptibility 



2„2 



Xt,q =r + 3|| + /3 q 



(CI) 



where < f3 < 1 and the same momentum cutoff A is 
taken for all three momentum components. This model 
naturally interpolates between the 2D case (/3 = 0) and 
the isotropic 3D case (f3 = 1). We show that the behavior 
of a c \ and a C 2 as function of /3 is very similar, although 
not identical, to the behavior of the two tri-critical points 
as functions on r\ z displayed in Fig. 13. 



The self-consistent equations for (j> and r become: 



r + u 



f A " dq K f 

Jo ^ % lo 



dq z X 



1 



+ 



1 



\r + ql+fPql-^ r + qf { + pq* + <t> / 
f A W dqu [ A * 

= § L Wo dqzX 



i 



i 



r + q 2 + I3 2 q 2 z - <fi r + q* + $>q% + , 



(C2) 



For simplicity, we set g — 1 below, but the results 
can be easily generalized for arbitrary g after rescaling 
02 1 jg_ We first do the 2D integral over g|| and then 
evaluate the one-dimensional integral in the q z direction. 
We define 



1(c) 



f Adq K f 

= Jo ^ qil Jo 



dq 2 



1 



ql + I3 2 q 2 + C 



(C3) 



with c = r ± <j). Evaluating the momentum integrals we 
obtain 

-^arctanf^ 



+Alog (A 2 + /3 2 A 2 + c) 
-A log (/3 2 A 2 + c) 



(C4) 



As before, we assume that A is large compared to both 
r and 4>. Then the first line does not depend on c and 
can be absorbed into the renormalizion of ro, i.e. into fo. 
The 2D and 3D results are indeed reproduced: in the 2D 
limit (3 — > 0, and we reproduce the logarithmic behavior 
and Eqs. (22); in the 3D limit (3 — > 1 and A — > oo such 
that arctan — > it/ 2 and we reproduce the characteristic 
square root behavior as well as Eqs. 27. 

A simple analysis shows that the crossover from 2D 
to 3D behavior occurs at a rather small /3 2 ~ 1/A. If 
a c \ and a c2 cross, they must cross in this regime. We 
rescale j3 2 by A (/3 — -^j=) and also rescale r and <f> by 

A (r = Af, <p = A(p), obtaining a cut-off independent 
equation in terms of (3 



4tt7(c) 2Vc 

— - — = — arctan 

A [3 




lof 



(/3 2 



+ 2 



(C5) 



The integral can be plugged into Eqs.(C2) to obtain the 
self-consistent equations for f and <f>. In order to obtain 
a C 2 corresponding to the nematic tri-critical point, we 
expand the self-consistent equations to second order in 
</> and look for the value of a when the coefficient of the 
</> 2 -term changes sign. A straightforward analysis yields: 
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a C 2 = 





/3V^ (3/3 2 + 5fea) + 


3( 




■ f c2 ) arctan 






P 2 r c2 (6 - 3/3 2 - 5f c2 ) + 3 (/3 2 + f c2 ) arctan (^L) 






(4 - /3 2 - r c2 ) 


+ 2 (/3 2 + f c2 J arctan (tJ=) 



(C6) 




as /3 increases, both a c i and a c2 increase, and their ra- 
tio approaches 1 from below. The difference between a c \ 
and a C 2, however, stays finite. 

We also solved numerically the self-consistent equa- 
tions (C2) for r and <j) and obtained the phase diagram 
in the variables fo and a for various /3. We found the 
same four phase diagrams as in Fig. 12. Namely, for 
small (3 we recover the behavior of Fig. 12 (a), while for 
j3 3> 1 (but with j3 < 1), we recover the behavior of Fig. 
12 (d). The two phase diagrams at intermediate /3 are 
also the same as those in Figs. 12 (b), (c). 



Figure 26: (Color online) The behavior of the magnetic and 
nematic tri-critical points q c i (dashed red/light gray line) and 
OLci (solid blue/dark gray line) as functions of the rescaled ft 
in the 3D anisotropic model with the bosonic susceptibility 
given by Eq. (48). Similar to the other anisotropic 3D model, 
Qfci and a C 2 cross at /3 w 1.7, and at larger /3, a C 2 becomes 
smaller than a c \. 



where f c2 is the critical value for the onset of Ising- 
nematic order: 



2 arctan 



1 = 



U0 



PV^2 



(C7) 



To calculate a c \ corresponding to the magnetic tri- 
critical point, we extend the self-consistent equations to 
include M in the standard manner and set r = <f>. The 
equation for the critical (\> c \ at which M = is 



2 v 7 ^ arctan ( 

'2<t>ci 



-log (^ 2 ) +log (/3 2 + 20 c 



(C8) 



Expanding the self-consistent equations in powers of 
M, it is straightforward to obtain the value of a — a c i for 
which the coefficient of the quadratic term M 2 changes 



sign: 



del 



arctan 



(C9) 



We plot a c i and a c2 as function of f3 in Fig. 26. We see 
that the two tri-critical points cross at a certain /3, be- 
yond which a C 2 becomes smaller than a c ±. As expected, 



Appendix D: Linear relationship between the 
nematic and magnetic order parameters 

In this Appendix, we show how near a simultaneous 
first-order magneto-structural transition, the nematic 
and magnetic order parameters obey an approximately 
linear relationship. We start with the self-consistent 
equations (41) in the magnetically ordered phase for ar- 
bitrary dimension 2 < d < 3. A straightforward manip- 
ulation leads to the equivalent equations: 



''0 



■4>(a-l) 
gM 2 



(Dl) 



with <j) = r. The condition dfo/dcf> — gives the value 
of Acf> for which fo is maximum and, therefore, the first 
instability of the system. Substitution of A(f> in the sec- 
ond equation then gives AM. Evaluating the algebraic 
equations yields: 



A<j> = 



(ff/2)*=* 



(d-2) 



AM = 2(g/2) 1 



a~ 1 
a (d-2) 



(D2) 



a — 1 



( 4 - d > /l-a(3-d) 
a- 1 



These are the values of the nematic and magnetic 
jumps for a < a c i — 1/ (d — 3). Notice that AM — > 
for a — > a c i and that the d — 3 result of Eq. (34) is 
reobtained from the second equation. 

Using the second equation of (Dl), we can expand it 
for (f> close to A(j> and M close to AM, obtaining: 

<j) - Acj) 4 [1 - a (3 - d)] 



A0 



(d-2) (a + 1) 
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M — AM\ a 2 /M-AMV 
AM J + V AM J 



+ {a-l)0 



M - AM \ 3 
AM ) 



with: 



2(d-2) (a + iy 



(D3) 



a 2 _ (6 - rf) (1 - 2a (3 - rf)) + a 2 (22 - 13d + 2d 2 ) 
ai 



(D4) 



Thus, the linear relationship between </> and M is al- 
ways present for small enough deviations from the jump. 
For a — > 1, this linear relationship dominates and extends 
to larger values of N. Indeed, as a — > 1, the jumps A</> 
and AM of Eq. (D2) become larger and, consequently, 

( "am^ ) becomes smaller for a fixed M. Furthermore, 

all the coefficients of the series expansion (D3) of order 
higher than quadratic go to zero, and the ratio ai/a 2 
(D4) between the linear and the quadratic coefficients 
tends to 2 for any dimension. 
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